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NAVAL ARCHITECTURE 


On long-term prediction of stresses in principal 
members of ship hull structures 


Current report 


Marian Bogdaniuk 
Polish Register of Shipping 


ABSTRACT 


The paper presents a discussion on the strength criteria and load cases required in the 
rules of classification societies for direct FEM strength analysis of ship hull structures, 
which show some differences to each other. The conclusion of it is that detail studies on 
stress values in ship structures are necessary to improve the requirements. With this end in 
view an effective method for long-term prediction of stresses in ship hull structure princi- 
pal members is proposed. The method is based on the concept of influence coefficients and 
spectral analysis of wave loads. A FEM model of the principal member system in the form 


of sD framed is applied to calculate values of the influence coefficients. Next, the concept is used of correla- 
tion factors for combining characteristic long-term stress values caused by global and local loads, with 
combined stresses due to general bending of the ship, zone bending of the principal members and local 
bending of longitudinals. As an example, results of stress prediction and correlation factors calculated in 
some points of hull structure of a panamax bulk carrier are presented and discussed. 


Key words : Rules for classification of ships, spectral analysis of wave loads on ships, dynamic stresses in ship hull structures 


INTRODUCTION 


Ship hull structures are usually designed according to safe- 
ty standards given by classification societies in their rules for 
the classification and building of ships, e.g. [1] to [6]. An im- 
portant issue of the rules are strength standards for the structu- 
res. Ship strength appraisal on the basis of such rules is perfor- 
med according to the scheme given in Fig.1. 


Calculation 
model 


Criteria 
(allowable 
stresses) 


Are the criteria 
fulfilled ? 
Acceptance 

of the structure 


Fig.1. Ship strength appraisal according to the rules 


Correction 
of dimensions 
of structural elements 


There are considerable differences in the rules [1] to [6] 
concerning design loads, and only in very few cases the [ACS 
member societies apply common requirements. Such case is 
the wave bending moment and shear force for hull general ben- 
ding in vertical plane, defined in [7]. 


An especially important issue for ship safety is strength of 
hull structure longitudinal members such as double bottom gir- 
ders, bottom and inner bottom stiffeners (longitudinals), etc. 

In these members normal stresses due to general bending 
of the ship, bending of the principal members (zone strength) 
and local bending of the stiffeners are superimposed (Fig.2). 


A-A 


web of girder 


E eater pias © neutral axis 
N of girder 


normal stresses 
due to zone 
bending 


normal stresses 
due to general 


deck bending 


bulkhead 


neutral axis 
of hull 


transverse 
principal 


web of girder 


bending longitudinal 

moment stiffener 
in stiffener of bottom 
due to local 

bending 


Fig.2. Normal stresses in hull structure longitudinal members 


The rules [1] to [6] require that strength of the principal 
members is to be appraised by applying FEM model of ship 
hull structure module at the midship region. The FEM model 
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should usually correspond to 3 successive holds. Some of the 
standard combined load cases are to be applied in the calcula- 
tions. The combined load cases are composed of global load 
components (bending moments and shear forces at general 
bending of the ship, for example) and local load components 
(external water pressure, pressure of the cargo, etc.). 

To compose the load cases, characteristic values of dyna- 
mic loads related to sailing in ocean waves, are multiplied by 
numbers from the range {-1;1}, called the correlation factors 
for combined load cases. The characteristic values of the loads 
are extreme values, i.e. usually those with the probability of 
exceedance equal to 10: 

For example, one of the combined load cases for strength 
appraisal of primary members, according to the requirements 
of DNV rules ([2 , 3]) and PRS rules ([4 , 5]), is to be compo- 
sed of the dynamic loads shown in Fig.3. 

heavy bulk cargo 
=) 0.59 M, 


empty empty 
0.59 My hold hold 


Ft ts f osp 


Fig.3. Dynamic load components acc. to DNV and PRS rules 


So, the correlation factor for the wave bending moment 
Mw takes the value 0.59, for the dynamic external water pres- 
sure - the value 0.5 and for the dynamic load corresponding to 
vertical accelerations ay - the value 0.5. 

Quite different values of correlation factors for the load case 
shown in Fig.3 are required by the rules [1] and [6]. According 
to [1] the value 1.0 is to be applied for My, , 0.4 - for ay and 0.5 
- for p. Whereas according to [6] the values : 0.625 for Mw , 0.0 
for ay and 0.5 for p are to be applied. 

Characteristic load values in the above mentioned rules of 
classification societies are usually defined at different probabi- 
lity levels and are to be calculated by applying quite different 
parametric formulae. Values of the allowable stresses are also 
different. 

These remarks allow to state that the determining of requi- 
red scantlings of ship hull structure members is based on very 
simplified assumptions. FEM analyses required by the rules 
would give quite accurate values of the stresses if values of the 
design loads corresponded to real loads accurately. The above 
given short comparison of some rule requirements of the clas- 
sification societies clearly shows that the required simply com- 
bined load cases must differ from the real loads. 

So, a.o. in the author's opinion the searching for more ac- 
curate design loads is the most important and urgent task for 
development process of the rules. Therefore in PRS a research 
task has been undertaken to find long-term stress distributions 
in longitudinal members of typical ship hull structures. Results 
of long-term stress predictions have been compared with stress 
values calculated as response of the structure to the design lo- 
ads defined in [4]. Analysis of the obtained results will make it 
possible to modify and improve the requirements concerning 
design loads [4]. 


LONG -TERM PREDICTION OF STRESSES 


IN LONGITUDINAL MEMBERS 
OF SHIP HULL STRUCTURES 


For many years PRS has been developing theoretical mo- 
dels and software for prediction of ship's motions and loads on 
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ship sailing in waves. The calculations based on the linear the- 
ory are performed according to the scheme shown in Fig.4. 


Calculation of transfer 
functions of ship’s 


Short-term 

prediction 

for motions 
and loads on ship 


Long-term 
prediction 
for motions 
and loads on ship 


response amplitudes 
to excitation from 
regular waves 


Fig. 4. Block diagram of the method for prediction 
of motions and loads on ships 

The transfer functions are calculated by using WAVE3D 
computer program developed in PRS. Three-dimensional wa- 
ter flow around the ship is taken into account. 

Short-term prediction calculations take the form of spec- 
tral analysis in which Pierson-Moskowitz spectrum in terms of 
the significant wave height and the mean wave period, is ap- 
plied [11]. 

The spectrum is narrow-banded and probability density 
function for the wave maximum values (peak values) takes the 
form of Rayleigh distribution, [8, 9, 10]. 

Long-term prediction concerns the whole life span of a ship 
(usually 20 years). The results of short-term predictions and 
statistical data in the form of joint frequency of significant wave 
height and mean wave period are used. The data for North Atlan- 
tic are used as usual. Detail description of PRS procedure for 
long-term prediction is given in [8] and [9]. 

The rules of classification societies should allow ship desi- 
gners to assess characteristic values of stresses in the hull struc- 
ture as effectively and accurately as possible. The characteri- 
stic values are those exceeded with sufficiently small probabi- 
lity value - e.g. 10. 

So, before formulating design loads and combined load 
cases for rule requirements one should start with analysis of 
long-term stress distributions in hull structures. 

The problem is that stress values in a transverse cross-sec- 
tion of primary member depend significantly on the loads ac- 
ting on a rather large portion of the hull structure because the 
primary members are mutually connected. 

To make the long-term prediction of stresses in primary 
members as effective as possible the influence coefficients of 
the loads on stress values at selected points of the primary 
members are to be calculated before the calculations according 
to the scheme shown in Fig.4 start. The concept of influence 
coefficients is explained in Fig.5. 


a? 


i-th piece 
of plating 


j-th piece 
of surface 
loaded by cargo 


Fig. 

The surface of ship hull plating loaded by the external wa- 
ter pressure p is divided into Nọ relatively small pieces under 
the assumption that p = const on each piece. 

Similarly the areas of inner bottom and decks loaded by 
pressure of cargo are divided into Ng pieces. Now the stress 
value o at a point of principal member can be calculated from 
the formula : 


5. The concept of influence coefficients 


Paj 


N N 
p a m 
i=l j=l 


(1) 


where : 
Np» Na - as defined above 


Pi - external water pressure at i-th area 

avj - vertical acceleration at j-th area 

Aj - the influence coefficient of pi on o 

Bj - the influence coefficient of Cj - ayj on © 

Cj - a coefficient to calculate paj (see Fig.5) 
corresponding to ayj (cargo density multiplied 
by height of the cargo). 


The values of A; and Bj can be calculated by applying 
a FEM model of hull structure and assuming p; = | or paj = 1 
on individual pieces of external plating surface, inner bottom 
or deck area, and zero values at the other pieces. The FEM 
model in the form of a 3D frame is sufficient for the system of 
primary members. 

To calculate transfer functions for amplitudes of o, the am- 
plitudes of p; and ayj (the complex numbers) calculated for the 
unit amplitude wave should be put into (1). 

The next steps of the calculations are to be performed ac- 
cording to the scheme shown in Fig.4. In the double bottom 
girders some components of 6 are superimposed (Fig.2). 

The phase angles of individual components of o for the 
ship on regular wave are usually different. This means that the 
characteristic values of the components can not be added di- 
rectly. The characteristic value of o corresponding to two com- 
ponents with characteristic values ©] and 62 can be written in 
the following form: 


o = Max(0,,0,)+c-Min(o,,0,) (2) 
where : 


Max (01 , 62) - the greater value of 0; and 62 
Min (61 , 62) - the smaller value of c1 and o2 
c -the correlation factor for combination 

of 61 and 62 ; this is a number from 

the range {-1 ; 1}. 


Such form of (2) corresponds to the above described me- 
thod of combining global and local loads according to the rule 
requirements of classification societies, where the correlation 
factors are used. 


RESULTS OF EXAMPLE CALCULATIONS 


Long-term prediction of normal stresses © was performed 
at 8 points of double bottom girders of a panamax bulk carrier, 
schematically shown in Fig.6. 


centre of short hold centre of long hold 


stool 
of bulkhead 


Fig. 6. Points in bulk carrier structure where long-term 
prediction of normal stresses o was performed 


Points 1 to 4 belong to the side girder 5.6 m distant from 
the ship's plane of symmetry. Points 5 to 8 belong to the central 
girder. 

The influence coefficients of external water pressure and 
cargo inertia forces on © were calculated by applying (1), and 
using the FEM model for the hull module principal members 


between the middles of two successive holds (see Fig.6), in the 
form of 3D frame. 

Symmetry of frame displacements was assumed at the ends 
of the model, i.e. girder's transverse cross-section angles of 
rotation around the axis perpendicular to ship's plane of sym- 
metry were assumed equal to zero. 

Full draught of the ship d=12.2 m was assumed. The cargo 
inertia forces corresponding to static cargo pressure on the in- 
ner bottom equal to 100 kPa were assumed constant over the 
whole inner bottom area. The long-term prediction calculations 
were performed by applying the above described method for 
the following locations of the FEM model along the ship : 


A -the transverse bulkhead located at the midship 

B -the transverse bulkhead 50 m distant from the midship, 
in the bow part of the ship 

C - the transverse bulkhead 50 m distant from the midship, 
in the aft part of the ship. 


The c-factor values calculated acc. to (2) for values of ©, 
6, and 62 at the probability of exceedance = 10°, at points | to 
8 (Fig.6) are listed in the table. 


Calculated values of c-factors 


C a | a f a O 
paps [ei [ee [me [ole 


The symbol A-1 in the table means the first case of 6; and 62 
combination for the location A of the FEM model, and so on. 


For the locations A, B and C of the FEM model the follo- 
wing combinations of normal stresses were considered : 


1. O0L=0p , 02= Oa 


where : 

Op - the stresses in the girders, at the level of the bottom, 
due to zone bending by external water dynamic pres- 
sure 

Oa - as above, but caused by cargo inertia pressure. 


2. 01 = Opa > 02 = OM 


where : 

Opa- the stresses in the girders, at the level of the bottom, 
due to zone bending by external water dynamic pres- 
sure and cargo inertia pressure acting together 

Om - the stresses due to general bending of the ship in ver- 
tical plane. 


3 0; =O, , 02= 0; 


where : 
Op - as at p.1, but at the level of bottom stiffener flange 
©; - the local bending stresses in the flange of bottom stif- 
fener, where : 


1? 
©] = ps 
12W 
s,l - see Fig. 2 


W - section modulus of the stiffener 
with strip of the plating. 
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CONCLUSIONS AND FINAL REMARKS 


Values of c-factors listed in the table suggest that combina- 
tion of global and local stresses o in hull structure longitudinal 
members is a complicated problem. 


The following interesting features 
of such combination method can be observed : 


a. Thec-values for double bottom side and central girders with 
the same x-coordinate along the ship are almost the same, 
in general. This suggests that the method of stress combi- 
nation in the form of (2) is reasonable. 

b. There are considerable differences between c-values for 
individual cases of the combination (A-1, A-2, A-3, B-1, 
etc.). It is interesting that the values of |c| in the cases A-1 
and C-1 are rather small. 

c. The c-values considerably depend on x-coordinate of the 
points in the girders. 

d. The c-values can be positive or negative. This is logical 
because the combination of stresses was considered. If com- 
binations of global and local loads which cause, at a point 
of hull structure, stress © equal to the result of long-term 
prediction based on direct calculations, were considered, 
then correlation factors would be the numbers of the same 
sign — independently of a position of a considered point. 


O The paper deals with the problem of combining the stres- 
ses. According to the rules of classification societies global 
and local loads are to be considered to create the load cases 
for FEM calculations. The above listed features b. and c. of 
stress combination suggest that the load cases required in 
[1] to [6] created by means of rather simple algorithms, can 
lead to prediction of stress values considerably different 
from the real ones. 

O The designing of hull structures on the basis of the rule 
requirements means that simplified loads are applied. The 
FEM calculations give quite accurate results for these 
loads, but calculated values of stresses can differ conside- 
rably from their real values. 

O For this reason in PRS a research task aimed at formulating 
the load cases to be used in the rules [4] in order to obtain 
more accurate results of stress prediction in longitudinal 
hull members has been carried out. The first step of the 
task is to gather a set of information on stress combina- 
tions. Then the load cases which cause sufficiently accura- 
te stress values in the structure will be searched for. 
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NOMENCLATURE 


ay - vertical acceleration 

A; , Bj - influence coefficients of pressure on normal stress 

c - correlation factor for normal stresses 

Cj - coefficient to calculate pressure values from acceleration 
values 

- spacing of floors 

- wave bending moment 

- number of pieces to which the inner bottom area is divided 

- number of pieces to which the area of plating is divided 

- external water pressure 

- spacing of stiffeners (longitudinals) 

- section modulus of stiffener cross-section 

- normal stress. 
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~ C onference = 


Scientific Meeting 
of 
Regional Group 


On 25 March 2004 the Mechanical and Electric Engi- 
neering Faculty of Polish Naval University, Gdynia, hos- 
ted the first-in-the year scientific meeting of the Regional 
Group of the Section on Exploitation Foundations, Ma- 
chine Building Committee, Polish Academy of Sciences. 


During the seminar three papers were presented : 


* Recycling problems in shipbuilding - research state- 
-of-the-art, issues to be solved — by W. Jurczak (Polish 
Naval University) 


* Research on influence of fuel charge nonuniformity of 
combustion engine on spectrum of torsional vibration 
of shafting — by St. Bruski(Polish Naval University) 


* Research on compression process to be applied to dia- 
gnosing ship piston combustion engines — by M. Lu- 
towicz (Polish Naval University) 


After interesting discussion the participants were acqu- 
ainted with a modern ship navigation bridge simulator and 
computer stands serving as training labotratories for navi- 


ae officers. A 
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On hydrodynamic forces acting 
on the ship in large motions 


Witold Błocki 
Gdańsk University of Technology 


ABSTRACT 


been applied. 


Present state of mathematical description of ship dynamic non-linear behaviour is pre- 
sented in this paper with a view to avoiding excessive complications in solving the pro- 
blem. The non-linearity concerns first of all Froude-Krilov forces and damping forces 
occurring after entering ship's deck into water or those resulting from drag of bilge keels. 
And, to the remaining, accompanying and diffraction forces the linear extrapolation has 


Key words : hydrodynamic forces, ship motions of large amplitudes, added forces 


INTRODUCTION 


The paper presents current state of mathematical descrip- 
tion of ship dynamic non-linear behaviour. The subject matter 
concerns large-amplitude motions of ship in waves, conside- 
red within the frame of ship roll non-linear theory. The split 
into the accompanying forces, Froude-Krilov forces and dif- 
fraction forces is assumed still valid. The Froude-Krilov for- 
ces are calculated with taking into account changeable wetted 
area of ship hull surface. However to calculate the accompany- 
ing and diffraction forces a method of extrapolation of linear 
solutions is applied. 


FROUDE-KRILOV FORCES 


Froude-Krilov forces are obtained by integrating pressures 
in waving water not disturbed by ship's presence in it. In the 
non-linear theory the integration is performed over wetted sur- 
face area which results from an instantaneous position of ship 
hull relative to wave surface. The pressure p is defined by the 
following expression : 


P-P, =pg(z—Ce*”) (1) 


The Froude-Kriloy resultant force F and moment M 
can be obtained with the use of the integration formulae : 


p-f- -7,7 p 
F= JEG p,)dS (FEHZZIav 


The forces are usually calculated by means of the above pre- 
sented surface integrals. In practice the wetted area is divided 
into a finite number of the directed area elements n AS, and 
integration is replaced by summation. The integration is rather 
troublesome as to determine the normal versor for every ele- 
ment AS is necessary. However the sufrace integrals can be re- 
placed by scalar volume integrals in accordance with Gauss- 
-Ostrogradski formula, which make determination of the normal 
versors not necessary. A way of calculation of the Froude-Kri- 
lov volume integrals was presented in [1], and realized in [2]. 


ADDED FORCES 


The added forces acting on ship result from the forced oscil- 
latory motion of ship in still water. They are determined by 
means of the formula : 


R=-m,ii;-Nji; ij=12..6 ® 


In the formula the Einstein's convention of summation is 
applied. F; , fori = 1, 2, ... 3, stands for x, y, z components of 
added forces, and for i = 4, 5, 6 it stands for components of 
moments of added forces (relative to ship's centre of gravity in 
the ship-fixed coordinate system). For j = 1, 2, 3 the velocities 
ùj stand for the components vx , Vy , Vz of the oscillation velo- 
city of ship's centre of gravity relative to its mean location, 
whereas for j = 4, 5, 6 the velocities stand for the ship's angular 
velocity components Wx , Wy , Wz (defined in the moving sys- 
tem). The accelerations tj for j = 1, 2, 3 stand for the accelera- 
tion components ax , ay , az of ship's gravity centre oscillations 
relative to its mean location, whereas for j = = 4, 5, 6 they 
stand for the components ©x , Oy ss 

The added masses m;j and damping coefficients Njj, called 
the hydrodynamic coefficients, are constant provided the ship 
in question performs oscillations of constant frequency. Then 
they functionally depend on the oscillation frequency w. If ship's 
motion is non-harmonic the added masses and damping coef- 
ficients do not have any constants and using them is senseless. 
In such case the hydrodynamic coefficients can be determined 
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by using the Fourier transform method. To this end it is enough 
to consider the addede force component F; resulting from j-th 
DOF oscillation as the Fourier transform is a linear operation. 


Assuming that the ship hull oscillations 
occur at the amplitude u4 : 


u(t) =u, sinot (5) 
one obtains the added force : 
F(t) = uo? m(@)sin@t—u,@N(@)cos@t (6) 


The hydrodynamic force can be expressed by means 
of the function of response to the step excitation [3] : 


F(t) = -fr (t)u(t — t)dt -fr (mült-r)dt O) 
0 0 


The relationship is also valid for harmonic motion. It is 
possible to introduce (5) to it, and next to compare it with (6) : 
o0 ice) 


m(@) = fi (t)cos@tdt — ea fr (t)sin wtdt 
OO 


0 
” z (8) 
N(@) =@ fo (t)sinœotdt + fr (t)cos@tdt 
0 0 
Moreover, the following relationships are valid : 
limm(®)=m,, lim N(@) = 0 (9) 
orn 


oon 


The function rz (T) must have the following forms in order 
the relationships (8) to have the limits given by (9) : 


r(t) = m,,8(7) 


where : ò (T) - Dirac's delta functions. 


(10) 


However the function rj(t) must be limited and decaying 
along with T increasing. Then the integral : 


co 


fr (T) cos wtdt 

0 
decreases along with  increasing,with the rate of about 1/o 
(its exact decaying rate is unknown). On accounting for the 
above mentioned comments and introducing the notion : r (T) = 
= rı(T), the formulae (8) are transformed into the following : 


m(®)=m,,— = r(t)sin@tdt 
0) 
e (11) 
N(®) = fa cos@tdt 
0 


The inverse Fourier transform makes it possible 
to achieve an unknown response function r(T) : 


r(t)= 2 flim- m(o)]osin otdo 
a 


2? (12) 
r(t=— [No cosætdo 
T 
0 
From the second formula (12) one obtains 
the relationship for the initial value of r (T) : 
2 o0 
r@=0)=2 [N(o)do (13) 
T 
0 
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From (12) it results that knowledge of one of the hydrome- 
chanical functions : m(@) or N(@) within the whole frequency 
range is sufficient to find the response functions r (T). Also, it 
can be observed that after determination of one of the hydro- 
mechanical functions and calculation ofr (T) the other function 
can be determined by means of the relevant formula (11). 


It is also possible to determine radiation forces in harmonic 
motion. On the basis of (7) the following is yielded : 
wo 
FO=-m ült- [r(a(t—-dt a4 
0 

The formula (14) is valid for any DOF. The integral term 
appearing in the formula is the so called memory effect which 
fast decays along with T increasing. Its decaying rate depends 
on that of r (T), hence it depends on the decaying rate of the 
damping coefficient N(@) along with œ increasing. The above 
presented theory, called Cummins model, is generally related 
to irregular waves. 

The formula (14) is generally valid for small harmonic 
motions of ship, performed around its equilibrium position in 
still water, within linear range. The linearity assumption is prac- 
tically valid until the bilge does not emerge from and the deck 
does not immerse into water. This makes it possible to extend 
validity of the formula into the case of large motions, as for 
practical reasons it is sufficient to be limited only to large roll 
amplitudes and small amplitudes of the remaining motions. 
Different calculation procedures for added forces at large he- 
eling angles can be met (Fig.1). 


A 


Fig. 1. Ship motions possible to be performed 
around its equilibrium position in still water 
The added masses mj; and the damping coefficients Nj, 
are numerically determined for the different draughts z and 
heeling angles ® [4] : 
mj; = mj;(Z,®) 


Na = Ni, (z,D) 


The forced oscillatory motion shown in Fig. la is more fa- 
vourable because the heeling angle ® is of a negligible influ- 
ence on added forces. 


(15) 


WATER ON THE DECK 


It is assumed that the added forces determined by (14) can 
be also applied to heeling angles greater than the angle of deck 
entrance into water. However in this case the additional forces 
due to presence of water on the deck should be taken into 
account. Unfortunately in the existing computer programs the 
problem has been neglected though results of model tests have 
indicated that the additional forces due to water on the deck 
would be significant [7]. It is only possible to roughly estimate 
them as the problem is insolvable strictly. 

The forces in question clearly depend on direction of mo- 
tion. If the deck emerges into water they practically equal zero, 
and if the water flows out from the deck they become signifi- 
cant. The forces can be estimated by applying the momentum 
conservation law, and the way of their estimation is highligh- 
ted in Fig.2 [8]. 


Fig. 2. Mutual relations between the deck area element dA, 
the normal to the deck, Ì, and the deck area element velocity V 


The momentum conservation law was used along the di- 
rection n for the water head over the deck area element dA. 
The water is subject to a change of momentum with time, which, 
according to 2nd principle of dynamics, is equal to the deck 
reaction Ra : 


d 
£L (mv,)=R, (16) 
rr n) 
or: 
ony, move =R, (17) 
dt dt 


The continuity equation for the water head over the deck 
area element dA subject to momentum change, can be written 
as follows : d 

ae pv, dA 
dt 
Taking into account the following 


auxiliary relationships (Fig. 2) : 


(18) 


v=b® , Un =V cos 


l . d : 
Vn =b Ọ cosð=l1ğ , 7 =] 


one can transform (17) into the following : 


R, =p®??dA+p@ldv 


The equation (19) expresses a unit force acting on the ship 
deck during water running off the deck. The elementary mo- 
ment of the force Ry relative to x - axis is expressed by the 
following relationship: 


dK = p°’ dA + pË’ dV 


(19) 


(20) 


The additional moment due to water running off the deck 
can be obtained by integrating (20) over the wetted deck area 
A and the volume of the water appearing on the deck, V : 


K =@’p[PdA + &pfPdv (21) 
A V 


The approximate formula (21) for the additional moment 
due to water running off the deck can be presented as the sum 
of two elements : the damping moment and the added mass 
moment : 


Oa g 
K =1,67 +i, @ (22) 
where : 
I, - 3" order moment of the wetted deck area 
i, - mass inertia moment of the water appearing on 
the deck, respective to the ship central axis x. 


DIFFRACTION FORCES 


The diffraction forces result from hull-induced disturban- 
ces of pressure distribution in waving water. In the linear the- 
ory the forces in question are calculated on the basis of wave 
action on motionless ship. The assumption is valid in the case 
of small motions only. Nevertheless it is also extended into 
large motions. The diffraction forces acting on a ship in irre- 
gular waves at large motion amplitudes are determined by 
superposing the forces resulting from particular harmonic 
components. It results from the assumption, common for ad- 
ded and diffraction forces, in which the ship is further consi- 
dered, despite large amplitudes of motion, as a linear object 
respective to those forces. The assumption makes calculations 
much simpler as it allows for using the characteristics of the 
ship in upright position, which are independent on time and 
instantaneous positions of the ship. Otherwise it would be 
necessary to determine them in every time step and for every 
instantaneous wetted area of hull surface, which is a very dif- 
ficult task. 

If to denote, by Fpsij and Fpgij, respectively the sinusoidal 
and cosinusoidal parts of the amplitude of the diffraction force 
associated with i-th DOF and resulting from j-th harmonic com- 
ponent, then the generalized diffraction force Fp; can be expres- 
sed as follows: 


Fpi(t)= >) 


j 


Fp ij COS(@pjt — €j) + as 
where the indicated summation is performed over all harmo- 
nic components. 


CONCLUSIONS 


¢ The radiation and diffraction forces are smaller than Fro- 
ude-Krilov ones approximately by one order of magnitude 
[5]; they influence first of all the phase shift angles betwe- 
en waves and ship motions, and — to a much smaller degree 
— amplitudes of the motions [6]. 

# On the basis of the subject-matter literature it can be stated 
that the differences between the linear theory of ship mo- 
tions and the non-linear one are not large, and that they do 
not significantly influence solutions of ship motion equa- 
tions if only motions of water relative to ship are assumed 
to occur within the range of ship's sides. 


NOMENCLATURE 


a - acceleration component of ship gravity centre oscilation 

A - ship deck wetted area 

b - distance between the deck surface area element dA and 

: the ship gravity centre G 

F - hydrodynamic Froude-Krilov resultant force acting on 
ship hull 

F; - components of added forces 

Fpi - diffraction force component 

Fpcij  - cosinusoidal component of diffraction force amplitude 

Fpsij  - sinusoidal component of diffraction force amplitude 


g - gravity acceleration 
i,j,k - versors 


k - wave number 

K - additional moment resulting from the water running off 
the deck 

l - distance between the deck surface area element dA and 
the ship plane of symmetry 

m - mass of water on the ship deck surface element dA 

mij - ship added masses 


- added mass corresponding to the angular frequency œ 
tending to infinity 
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Kd 


- hydrodynamic Froude-Krilov resultant moment acting on 
ship hull 

- unit vector normal to ship surface 

- ship damping coefficients 

- pressure inside wave 

- atmospheric pressure 

- response to the step excitation 

- tracing vector 

a - unit force acting on the deck surface area element dA 

- instantaneous ship hull surface wetted area closed by free 
surface of wave 


PO Zp 


D 


nAWrH 


t - time 

u - oscillation of ship gravity centre 

uA - ship hull oscillation amplitude 

uj - oscillation of ship gravity centre in the j-th direction 


respective to its mean position 
v - velocity of the deck surface area element dA 
Ui, - projection of the velocity V onto direction of the normal n 
V - instantaneous volume of immersed part of ship hull; 
volume of water on the deck 
x,y, Zz - spatial coordinates 


£j - phase shift angle of harmonic component 
E - wave profile ordinate 
0 - angle between vectors of the velocity D and the normal n 
p - water density 
T - time shift angle 
p - ship heeling angle 
ğ,ğ - ship roll angular velocity and acceleration, respectively 
(0) - oscillation frequency 
Op - encounter frequency of the ship 
Acronims 


CTO - Ship Design and Research Centre 
DOF - degrees of freedom 

ISC - International Shipbuilding Conference 
PRS - Polish Register of Shipping 


On 18 February 2004 the 10th, successive Seminar on 
Drives and Control Systems 


was held in Gdansk. The series of such seminars accompa- 
nying the Fairs of Manufacturers, Subcontractors and Pro- 
viders of Power Units and Control Systems, have been each 
year organized by Gdansk University of Technology since 
1995. This has been a good occasion for direct mutual con- 
tacts of representatives of scientific centres working at de- 
velopment of this engineering branch and the manufacturers 
and dealers of such technical devices, which next brings about 
novel solutions based on results of R&D projects. 


As usual, the Seminar was devoted 
to the problems of the following areas : 


drives and control systems of machines 

automation and dynamic behaviour of driving 

systems for shipbuilding and power industries 
e electronic devices applicable 

to driving and control systems 

automation of electric drives 

applications of hydraulic and pneumatic drives 

applications of control 

and signal processing methods. 


( ‘onference 


10th Jubilee Seminar 
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The problems were discussed during three sessions : 


A - Mechanical, hydraulic and pneumatic drives 
(12 papers) 

B - Automation of electric drives (6 papers) 

C - Applications of control and signal processing methods 
(7 papers), 


and, during one poster session 
another 22 topics were presented. 


Moreover, the Seminar was accommpanied with the me- 
etings organized in the frame of the workshops devoted to 
the following topics : 


A I -Control techniques for linear electro-hydraulical 
servodrives. (Technical University of Poznan) 

A II - Advanced simulation techniques for converter dri- 
ves. Application of Tcad package. 
(Gdansk University of Technology) 

à JII- Automation of ship electric power system. 
(Gdynia Maritime University) 

A JV- Electronic microsystems - their design, diagnostics 
and integration. (Gdansk University of Technology 
in cooperation with TASC company). J 
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From 1 to 4 June 2004 a group of members of the Ger- 
man Society of Naval Architects and Marine Engineers 
(Schiffbautechnische Gesselschaft - STG) took part in its 
yearly summer meeting held abroad. An important event of 
the meeting has always been a scientific conference. This 
year it was held on 2 and 3 June at Szczecin Maritime Uni- 
versity, the local host of the conference. 


The conference was carried out in four sessions : 


ve Inland Shipping; Strength — 3 papers 

ve Hydrodynamics — 4 papers 

ve Ship Design and Ship Safety — 3 papers 
ve Ship Propulsion — 8 papers. 


During the last session 4 papers were presented 
by scientific workers of Szczecin Maritime University : 


> An experimental investigation of influence of fuel injec- 
tion timing on exhaust emission of marine low-speed 
engines — by T. Borkowski, D.Sc. 

> Examinations of structure and selected chemical proper- 
ties of the sludge produced during cleaning process of 


( ‘onferences 


STG Conference in Szczecin 


residual fuel, lubrication oil and bilge water 
by P. Rajewski, D.Sc. 

> Design and facilities of a testing stand of oily ship water 
cleaning process — by J. Listewnik, Assoc. Prof., D.Sc., 
A. Wiewiora, D.Sc. 

> Identification of marine generators in service and main- 
tenance — by P. Bielawski, D.Sc. 


The Conference participants had the opportunity of visi- 
ting Szczecin including its historical monuments and a ship- 
yard, as well as to take part in tours to neihbouring small 
historical towns and health resorts. 


Moreover the host of the Conference 
presented some of the University's research facilities 
and research projects, namely: 


O Newly developed homogeniser 
and its performance tests 

O EU- funded research project 
on machinery maintenance 

O Bilge water deoiling test stand. 


A conference under sail 


From 28 March to 3 April 2004 Ist Domestic Scientific 
Conference, an original conference on : 


Scientific and technical problems 
of professional sailing 


was held on board the sail ship POGORIA 
during her voyage on the route : 
Genoa - Elba - Bonifacio - Nice - 
- Monaco - Portoferraio - Genoa. 


It was organized by Faculty of Mechanical Engineering of 
Gdynia Maritime University, Faculty of Mechanics, Energy 
and Aeronautics and Faculty of Motor cars and Heavy Ma- 
chinery of Warsaw University of Technology, and had on its 
agenda presentation of 12 papers split into 3 topical groups : 


Measuring and testing 


** On possible assessment of a yacht in natural conditions 
by Zb. Dabrowski (Warsaw University of Technology) 
(an invited paper) 

“* Multi-channel data acquisition and processing system 
for steering the sailing boat by means of TCP/IP 
by M. Pilarski (Warsaw University of Technology) 

* Design of a system for determination of polar characte- 
ristics of yacht speed — by P. Ktopotowski (Warsaw Uni- 
versity of Technology) 

% Application of MLS method in testing acoustic qualities 
of small compartments on yachts — by W. Batko, T. Wszo- 
tek, W. Wszotek (Mining - Metalurgical Academy) 

% Vibration propagation measure as an index of yacht mast 
profile weakening (abstract) — by G. Klekot (Warsaw 

C University of Technology). 


Novel materials and engineering processes 


* Some problems of manufacturing and operation of com- 
posite masts (an invited paper) — by W. Skorski (Warsaw 
University of Technology) 

Dynamic features of modern composite masts 

by P. Detuszkiewicz, J. Dziurdz (Warsaw University 

of Technology) 

Problems of modelling the polimer - structure materials 
by M. Dudziak 

Problems of transferring the concentrated forces into 
composite structures — by M. Rodzewicz (Warsaw Uni- 
versity of Technology) 

Problems of application of metal sandwich structures for 
elements of sea-going ships (abstract) — by J. Kozak 
(Gdansk University of Technology). 


* 


* 


Self - steering gear and autopilot 


Design of a self-steering gear having an actively control- 

led vane-fixing angle — by P. Gorecki, B. Stankiewicz, 

P. Szajkowski (Warsaw University of Technology) 

A simple yacht autopilot controlled by means of electro- 

nic compass — by P. Ktopotowski, St. Misiaszek 

(Warsaw University of Technology). 
7 = 7 =EN BE 
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Roll response of ship fitted 
with passive stabilizing tank 


Ludwik Balcer 
Gdańsk University of Technology 


ABSTRACT 


Physical and mathematical models of roll motions of a ship equipped with a roll stabili- 
zing tank of working liquid free surface, is presented. Elaboration of the physical model 
was based on the idea of two mutually coupled mathematical pendulae. On the basis of 
the physical model, motion equations of the ship with the tank were determined and sol- 
ved. A way of using the achieved solutions is shown, as well as calculation formulae for 
coefficients of the motion equations, directly related to the main parameters of the ship 
and tank, are presented. Such form of the coefficients enhances possiblility of application 


of the equations and their solutions in ship design practice. Some examples of the use of the solutions for 
analysis of stabilizing effectiveness of a designed tank for a given ship, are also attached. Moreover, guide- 
lines for correct design of the stabilizing tanks having free surface of liquid, based on the proposed physical 
model of the ship-tank system, are offerred. It is also indicated that on the basis of the presented results it 
would be possible to search for ways to make operation of the stabilizing tanks in question more effective. 


Key words : seaworthiness, ship safety, ship hydromechanics, stability of floating units 


INTRODUCTION 


Out of the oscillation motions performed by a ship in wa- 
ves the side roll (heeling) undergoes stabilization most often. 
The stabilization means first of all a limitation of roll amplitu- 
de values, sometimes associated with increasing its period. The 
stabilizing tanks of free surface of liquid, which are further 
considered in this paper, are simple devices intended for limi- 
ting roll motions. 

A disadvantage of the existing methods describing work of 
stabilizing tanks is their low usefulness for designing. The equa- 
tions applied in the methods are rather complicated [8, 10, 12, 
13 , 14] as they are usually aimed at correct describing real 
motions of ship fitted with tank, and not at achieving practical 
design recommendations. 

In this paper is presented such physical model of the ship 
with stabilizing tank under rolling and mathematical model 
based on it, as to obtain motion equations of a possibly simple 
form more useful for formulating important design recommen- 
dations. Coefficients of such equations are directly associated 
with ship and tank parameters. The obtained form of solutions 
makes it possible to get access to these parameters of ship and 
tank whose selection is crucial for ensuring expected stabili- 
zing effects. 

The simplifications introduced to obtain an appropriate form 
of the equations and their solutions do not impair to an evident 
degree the quality of the description of ship motions, based on 
them. 
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FREE - SURFACE STABILIZING TANKS 


The basic types of free-surface stabilizing tanks presented 
in Fig.1 , 2, and Fig.3 and 4 exemplify their possible location 
on ships. An appropriate adjustment of natural period of mo- 
tions of liquid contained in tank to that of rolling ship or of 
wave acting on ship, is of a dominant influence on stabilizing 
effect. The way of changing the motion period of liquid in such 
tanks consists in changing level of the liquid, which also leads 
to changing its amount. 


© N palisade 
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Fig. 1. Types of cubicoidal stabilizing tanks [16] 
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Rys. 3. Location of "Flume" tank onboard a dry cargo ship 


—_ 


ilizing tank 


Rys. 4. Stabilizing tank in the form of container 
installed on deck of a fishing ship [2] 

Behaviour of a ship with stabilizing tank can be limited to 
consideration of motions of the system of two degrees of fre- 
edom (DOF), namely : for ship — side rolling angles, for tank — 
-translations of centre of mass of the liquid contained in it. The 
liquid mass centre translation can be represented by the average 
slope angle of its surface relative to ship, against its initial posi- 
tion. Hence description of motions of the ship-tank system amo- 
unts to the known problem of motion of a two-DOF system, 
whose solutions are based on the Lagrange 2nd kind equations. 
However in order to form such equations a "physical" model of 
the ship-tank system, i.e. physical representation of the ship, of 
the tank and their mutual coupling should be first determined. 
A form of external excitation of motions of the system, due to 
wave action on ship, should be also given. 


PHYSICAL MODEL 
OF SHIP - TANK SYSTEM 


A rolling ship is usually represented by a physical pendu- 
lum of an inclination moment dependent on ship's metacentric 
height. The mathematical pendulum is the simplest and ligh- 
test out of possible physical ones of a given inclination mo- 
ment and period. The system of two mutually coupled pendu- 
lae is a known mechanical analogy of the two-DOF oscillating 
systems. It was assumed that the ship with stabilizing tank can 
be represented by the system of two mathematical pendulae 
appropriately selected and mutually coupled. The main oscil- 
lation axis of the pendulae corresponds to location of ship's 
rolling axis. The location of ship's rolling axis was assumed 
constant, known and not necessarily identical with that of ship's 
mass centre [1]. A schematic diagram of the physical model is 
shown in Fig.7, whereas Fig.5 and 6 highlight its association 
with the main parameters of ship and tank. The necessary mo- 
delling principles according to which the elements of the pen- 
dulae and their coupling parameters have been determined, are 
given by the formulae (1) below. The main coupling parameter 
of the pendulae is the distance between their rotation axes; i.e. 
the segment Og m’ in Fig.7. The point O% represents the loca- 
tion of the ship's rolling axis, Os. 


Fig. 6. Ship with stabilizing tank, heeled by the statical angle Q 
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Fig. 7. Physical model of a ship with stabilizing tank 


Modelling principles 


The system of pendulae shown in Fig.7 represents a physi- 
cal model of a ship with stabilizing tank having free surface of 
liquid. The model can be deemed correct if statical and dyna- 
mical response of the system to external excitation moment is 
the same as that for the real ship with tank. For the pendulum 
of ship and tank the following can be written : 


2 T? g 
l= o Ig=—2 2) 
h, Aa? 
2 T? 
pei g jo = 6) 
a r, i 4T? 
where : 


Tæ - period of natural oscillation of ship's mathematical 
pendulum, and of ship itself 

lp -length of ship mathematical pendulum 

Tsm - radius of ship mass inertia 

Tym - radius of mass inertia of liquid in tank 

rz -metacentric radius of tank 

g -acceleration of gravity a 

hg -initial metacentric height of ship without tank, GMs 

Tp -period of natural oscillation of tank mathematical 
pendulum, and of motion of liquid in the tank 

lọ -length of tank mathematical pendulum. 


The static condition can be expressed by the relationships 
(4) and (5). The expressions (4) concern the equality of the 
moments due to heeling for the ship itself, the tank itself and, 
respectively, for the pendulae being their physical models. By 
using them the value of the mass mø - for the ship mathemati- 
cal pendulum, and of the mass mg - for the tank pendulum, can 
be obtained : 
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Modelling principles 
2 2 
ToS T8 
4T 4AT? 
m, h m,r 
Mo = — m, =—— 
i I, 
m 
R'=R — 
Mo 
fe —Mo8 
Mo = Dhamsin 0t 
(1) 
m, h 
m,h, = Mmo lọ , Mo = i i 
® 
a (4) 
m, t, = Mo l > m= Ls 
lo 


where : 


ms - ship mass (without tank) 
mz - tank mass, i.e. mass of liquid 
contained in the tank. 


The total static condition has to account for similarity of 
the righting moments for the heeled ship with tank and for the 
system of pendulae inclined by the same angle. Rigidity of the 
ship with tank and of the pendulae can be expressed by using 
the quantities shown in Fig.5, Fig.6 and Fig.7, which provides 
the following equality : 


msh, -mR = molo — MpR' =(m,;+m,)h-—mzr, (5) 


where : 


- the distance between the tank metacentre and 
the metacentre of the ship's immersed volume 
increment AV 

R'=m'c' - the distance, in the physical model , which cor- 

responds to R of the ship 

h - initial metacentric height of the ship with tank 

without any correction for liquid free surface 
mfz = ixp - correction for free surface of liquid contained 
in tank. 


R=mc 


To the system of pendulae the same moment due to water 
in tank must be applied as that applied to the ship. Knowing 
that : msh, = m@l@ one obtains, acc.to (5) : mR = mR’. 
It corresponds to the moment of couple of forces. For the ship 
with tank these are the forces : p; = f, = m,g (Fig.6), and 


respectively for the system of pendulae Po = fo’ = mo (Fig.7). 
Hence in the physical model the value of R' amounts to : 


R=R=2 , R=z,+1,—-T-AT/2-p, 
Mo (6) 
Pe = Alwx /AV 
where : 
Pe -metacentric radius of displacement layer 


of AT in thickness [7] 
Alyx - increment of transverse inertia moment 
of waterplane area for a given AT. 


The dynamic condition first of all consists in the equality of 
the inertia moment of the ship and that of its pendulum, as well 
as the mass inertia moment of liquid in the tank and that of its 
pendulum. Moreover has to be fulfilled the equality of kinetic 
energy of the ship with tank and that of the system of pendulae 
when the same roll inducing moments are applied to them. 


Hence the following is obtained for the ship and tank : 


Meron = mol g m2 = file (7) 
All the elements in that expression have been already de- 
termined by the statical condition. By means of the expres- 
sions (2), (3) and (4) it is easy to check that the above given 
relationship is satisfied. Kinetic energy of the ship with tank 
will be equal to that of the system of pendulae, if the following 
relationship is satisfied : 
msn + m,Rê + mn = mol + mR’ + Melé 
where : 
Ro - the distance between tank's metacentre and 
ship's rolling axis, mO, 
Rọ - the distance between the suspension point of 
the tank mathematical pendulum and the oscil- 
lation axis of the system of pendulae, m'O'; . 


From (7) it results: mR; = mR’ , hence, after taking into 
account (4), one obtains the following expression for Rọ : 


l r 
=R, | or: R} =R, = 
0 0 r 0 0 r 

Z zZ 


Ro = zp + rz—T-—AT + ay 

The quantity aw appearing in (8) determines the distance of 
the rolling axis from the waterplane (Fig. 5) [1]. It concerns the 
draft T and the mass centre location coordinate zg for the ship 
with "frozen" liquid in the tank. 

General information on effects of operation of a tank selec- 
ted for a given ship can be achieved on the basis of the roll 
trasfer function of ship with tank. Therefore the roll is induced 
by the heeling moment resulting from the action of regular plane 
beam wave on the motionless ship. Wave frequencies to be 
selected should properly cover the whole range of encounter 
frequencies possible to occur in service of the given ship. The 
roll inducing moment due to regular plane wave is of the follo- 
wing form [3,11]: 


(8) 


M,=Dhamsina,t , Om=Kg Kræ, (9) 
where : 
Mı - ship roll inducing moment 
D - ship buoyance force 
h - ship metacentric height 


QÀ , m - amplitude and effective amplitude of wave slope 
angle, respectively 

w - wave frequency 

t - time 

KB , Kr - wave slope angle corrective coefficients dependent 
on B/À i T/À , where A stands for wave length. 


In the expression (9) the buoyance force D and metacentric 
height A concern the ship with tank. It may take also another 
form if influence of sway on motions of liquid in the tank is 
accounted for. To take into account the sway it was assumed 
that it influence solely motions of liquid in the tank [15], and is 
equivalent to the horizontal component of orbital motion of 
the ship in real waves. The radius of the motion is equal to 
a half of wave height associated with the amplitude of its effe- 
ctive slope angle, Œm. And, for the horizontal oscillations the 
following is valid : 


H 
Val esr, aa) 


therefore : 


» Q, (t) = Qm Sino; t 


y,(t) = V, = ao OA mCOSO, t (10) 
Oy, 
where : 

yx -ship sway 
Hm - wave height corresponding to the 

effective slope angle Om 
V} - horizontal transverse component 

of ship sway-induced velocity. 


In the physical model the accounting for sway is equiva- 
lent to the applying of the horizontal oscillations V}, to the 
pivoting axis of the tank pendulum (Fig.7). The so obtained 
motion equations of the tank mathematical pendulum, at ne- 
glecting the damping, are as follows : 


H 
+ =o = spo sin@, t 


(0) o (11) 


" 2 2 $ 
POPE OFA, SiIn@,t 


It can be observed that instead to take into account ship 
sway in the physical model it is possible to apply, to the tank 
pendulum, the excitation moment whose normalized form is 
given by the right-hand side of the relationship (11). The same 
value of the moment was used in compliance with the Watana- 
be's method, in [15]. 


MATHEMATICAL MODEL 


The mathematical model is equivalent to the set of motion 
equations of ship with stabilizing tank, achieved on the basis 
of the above presented physical model. A searched solution is 
the roll transfer function further used for analyzing the beha- 
viour of ship with tank. Solutions are searched for within the 
frame of linear approach. 


Derivation of motion equations 


The Lagrange 2nd kind equations are used 
for derivation of motion equations : 


q, 2) 


In the equation the angle ® - for the ship's pendulum and 
the angle @ - for the tank's pendulum corresponds , in the phy- 
sical model, to the quantities qj , namely qı and q2 , respective- 
ly [9, 12, 13]. 

Particular kinds of energy considered in the physical mo- 
del are the following : E - kinetic energy, U - potential energy 
and  - dissipation energy. In Fig.7, the horizontal line passing 
through the axis of the system of pendulae represents the plane 
of zero- value potential energy of the system. 
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The kinetic energy of the system of pendulae is of the form : 
2 ob? 2 ob? 
E= Mo CEJ +m (lọ —Ro ) au 
)@V, + 
Vi 
” 2 


The potential and dissipation energies are as follows : 


t Mo (l= Ro 1, bo E Mo d= Ro 
-2 (13) 
+ml,@V, +m, tm 


2 
U = [moglo - mpg (R -1p)] 2 + 
32 (14) 
+ mp Blo Po + mg lo 


a S O 
©=—k,@0°+—k (15) 
5° 70°? 


The physical model of tank and ship represents a two-DOF 
object, hence the expression (12) is the set of two equations. 
One of them is the equation of derivatives of relevant energies 
respective to ®, (q1), and another - respective to @, (q2). The 
wave-induced moment Mg; = Mg; = Mj, complying with the 
equation, is introduced to the right-hand side of the ship roll 
equation. As the passive tank is considered, Mg; = Mq2 = 0 is 
introduced to the right-hand side of the equation of motions of 
liquid in the tank. After calculation of energy differentials, re- 
placements and ordering the equations, one obtains : 


[mols + mg (Rô — lp) ] ®+ko Ò + 
+ [mo g lo- me g(R’— lọ)] ® + 
-me (Rå — lo) 1G + Mp Flop + 
—Mpo (Ro — lẹ) V= D ha msin at 
Mo 14 6+ koG+ Mo g lP — Mg (Rô — Ip) lọ Ë + 
+ Mo g lọ D + my lọ V. = 0 


The final form of the equations is obtained by introducing 
the value V} = gOm sino3t into (16), with making use of (10). 


(16) 


Therefore : 
[mo lő + m, (Rò — lọ)” ] Ë + ko + 
+ [mo g lo- My g (R’— lọ)] ® -me (Rô — lọ ) lọẸ + 
+ mo glo @ = [Dh+mog (Rò — lọ)] a msin ont 
Mo lọ P+ ky + Mo g lọ P- Mg (Ro — 1p) Lye + 
+ Mo Z lo =- Mg lọ Z A m Sin Mat (17) 


Under the made assumptions the expression (17) represents 
the full form of the roll motion equations of the ship fitted with 
passive stabilizing tank. The only simplification represents not 
taking into account the mutual influence of ship's oscillations 
and motions of liquid in tank, due to damping. If an influence 
of sway is neglected the motion equations take the following 
form : 


[mol + mg (Ri — 1,)?] ®©+ko@ + 
+ [mo g lo- Mg g (R’= 1p)]  — My (Rô — lọ ) lọ + 
+ mo glog=Dhamsin at (18) 
Mo lé +koG+ meglo E+ 
-m (Rô — 1p) lt my gly =0 


Each of the coefficients appearing in the motion equations 
should be clearly defined by ship and tank parameters. The 
equations (17) and (18) derived on the basis of the proposed 
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physical model, satisfy the postulate. Their normalized form is 
as follows : 
os s 2 oe 
D+ Lo ® + Wo O- BiG + 
+ 100° = Opi Am sin ant (19) 
Ọ + Ho G+ -bË 3 Oy D= 


== Wp Om SIN Wy, t — accounting for sway 


Ö + to ® + Mo? D- BiG + 
+ y1Oo Ọ = ine” Om SiN œt (20) 
@ + Ho $+ 9 p-bÖ+ 


+ Wo” O=0 — not accounting for sway 


Next, the expression appearing at ® in the equations (18) 
was marked K and the following notation was applied : 


[mols +m, (Rå - ly)*]=K , koK = poi 
[mo g lo- m, g (R= lp)] K = o61 

Me (Ri — lo) lK =B , meglo /K = y1 WH 

[Dh+m g Rà- l] K=o51 , DhK=om 

ky /mg lg =o >» Mp Flp/Mels = g/ly = ws 

Mg(Rb — lo) lp /mg1§ = (Rè — lọ)/lọ =b 


The usually available initial data deal separately with ship 
itself and tank itself. Then mols is given istead of K. Hence 
the coefficients in the equations (19) should be expressed by 
means of : 


(21) 


2 
Har =52 3 ai a > B-E 
Y > _ O% 5 2 K” = 
Ne a Ml See 3 Mp2 = Po Fr 
where 
uo =ke /me lg > og = mo g lo/mo lő = g/l% 


The quantities ug, wd, B, y, K', K", K", expressed respecti- 
vely by the physical model parameters and the ship and tank 
parameters, have the same form as that given in the expres- 
sions (23), which makes it possible to examine structure of 
each of the coefficients appearing in the motion equations. 


PHYSICAL MODEL SHIP 
_ ko _— ko _ ko 
mg Mo lb ĝi L + lxx m, n 
megl m, gh, gh, 
Op = So = 8 Wo = : B 
Molo lo lily Gn 
m,| r 
y=—** y=E4 8) 
Mo lo s 
m, (Ro-!,) 1 R,k 1 
f=. 2 0 to) “9 B=y oo _] — 
mp lp lọ v lot, ko 
eee oe a 
ly lo £ 


PHYSICAL MODEL SHIP 
Ri -1,)° 2 
K’ 1 Mo ( 0 o) K’=14p? Ko 
Mo ls 
K”= “9 D K” = = 
a o jeg | ed 
Mo ly 1z 
geaen K"= a+ 
i 
Tz 
l 
R) = Rọ | Ro = Zp + rz— T — AT + ay 
Tz 


AT 


h = net (THAT + pe ~ zy hs) 


1+€ 


Pe =1s (1-8/0) 


— 2? _ 
2+(S/ajé 


E=m,/m, , ko=@p/@o , M0 =A/M , Led = o/Oo 


2 


aw = 0.43 [(T + AT) — zg] + 0.1B 


for : AT/2 ~0 and pe ~ 0 it yields : 


RZ +1r,—T | 
paite 
1+6 


Solutions of the motion equations, 
transfer functions 


The general form of ship roll transfer functions 
obtained as a result of solving the equations 


Ro = R + ay 


(T - zp- hy) 


(19) and (20), is the following : 


The presented expressions which make it possi- 
ble to calculate transfer functions for the ship fitted 
with stabilizing tank, directly provide the following 
options : 


+ to follow changes appearing in solutions depen- 
ding on an applied variant of the equations 

+ to determine in which way changes of the main 
parameters of tank and ship influence the calcu- 
lated functions. 


In the equations based on ship and tank data 
axcess to variables and parameters is easy and di- 
rect. It is possible to directly change such parame- 
ters as: Zp - tank position, kọ - tuning factor of natu- 
ral frequencies of tank and ship, r, -metacentric ra- 
dius of tank , i.e. its type, € - i.e. amount of liquid in 
tank, Ds and h; - i.e. ship loading state. 


The assumed form of the data and auxiliary 
values, (23), including the variables : 


w0 = 0)/W@ - dimensionless excitation frequency, 
an independent variable 
Uoo = Le/W@- dimensionless damping coefficient of 
motions of liquid in tank, the para- 
meter for sucessive realizations of 
functions of a given type 
ko = W/W - relative natural frequency of oscilla- 
tions of liquid in tank, possible to be 
used as an intermediate parameter 
- relative mass of liquid in tank, possi- 
ble to be used as an intermediate pa- 
rameter; 


= mz/Ms 


all of them make it possible to use a single, simple 
set of initial data, regardless of simplification level 
of the equations in question. 


APPLICATION OF THE 
MATHEMATICAL MODEL 


The obtained roll transfer functions make it possible to fol- 
low influence of the main ship and tank parameters on ship roll 
amplitudes, and thus to ensure stabilization effectiveness of 
a given tank. To achieve courses of relevant functions the data 
ofthe ship and tank considered in [16] were selected. A correct 
course of the achieved functions serves also as a check of cor- 


rection of the formulae based on the presented physical and 


mathematical models. 


a 


© A’ +B? 
Ks 2 2 
m C+D 


(24) 


For (19) where sway is accounted for, one obtains : 


A = poo Ko @0 g 
B =k? (1 +y) -00 (1 -B kð 


C=-K' (1/K"— @”)(ko” — 0”) + 


+y + Bao (ko? + b 00^) + Loo Hoo Ko wo” i 


D=K' Ho ko Wo (1/K"— o) + Uoo Mo (ko — 0°) ° 


K,=1 


For (20), without accounting for sway, it is : 


A= Hoo ko Oo 5 B= ke ag 


C and D as for (20) , K,=K"= 


m +m, h 


tig lo 


The ship with stabilizing tank acc.to [16] : 


Ship parameters : 


length b.p. Lop = 266.00 m 
breadth B = 42.50 m 
draught T= 10.85 m 
mass displacement A = m; = 97800 t 
height of centre of gravity ZGs = 15.40 m 


initial metacentric height 


(without tank) GM, = h, = 3.576 m 
block coefficient ô= 0.797 
waterplane coefficient a = 0.831 


natural roll period 
(natural frequency) 
roll damping coefficient 


to = 17.10 s (w@ = 0.37 1/8) 
uo = 0.070 1/5. 


Tank parameters (see Fig.1) : 


length 
breadth 


1, = 13.60 m 
b, = 40.89 m 
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- design height hy = 4.79 m 
- natural frequency of motions 

of liquid in tank We = 0.36 1/s 
- level of liquid in tank h, = 2.20 m 
- mass of liquid in tank 

at p = 1.00 t/m? m, = 1224 t 
- tank bottom height over 

the plane of reference Zq = 20.00 m 


- tuning factor 
- relative mass 
of liquid in tank 


ko = We/O@ = 0.97 


E = m/m, = 0.0125 


The diagrams shown in Fig.8 and 9 illustrate the roll trans- 
fer functions of the ship with tank. The dimensionless dam- 
ping coefficient of motions of liquid in tank, uoo, is a parame- 
ter for particular curves. Fig.8 deals with the full set of equa- 
tions in which sway is accounted for, and Fig.9 - the full set but 
without taking into account those oscillations. The presented 
diagrams are usually used for choosing an optimum value of 
the damping coefficient Uo, and after that — for evaluating the 
maximum stabilizing effect [5, 12, 15]. A value of the dam- 
ping coefficient for tank is so selected as in the middle, "de- 
sign" range of the diagrams the difference between the roll 
amplitude of the non-stabilized ship, which corresponds to the 
value [go = %, and a selected value e.g. Ugo = 0.6; were the 
greatest. The value of this difference, determined for the maxi- 
mum amplitude of the non-stabilized ship, hence for the exci- 
tation resonant frequency, is the stabilizing effect searched for. 
A value of Ugo should be selected in such a way as in a range 
far from resonance e.g. outside the characteristic nodes of the 
transfer functions' diagram, an increment of amplitudes, usual- 
ly appearing there for stabilized ship, were the smallest. The 
diagrams also show that for ugo = 0, which is possible theore- 
tically only, such excitation frequency exists at which the sta- 
bilizing effect amounts to 100%. The stabilizing effect expres- 
sed in percent, is the ratio of the roll amplitude reduction and 
non-stabilized amplitude. 


i] 
1.00 


0 
0.20 
Fig. 8. Diagrams of roll transfer functions with taking into account sway 


0.40 060 0.80 1.20 1.40 1.60 
Fig.8 it can be observed that when sway is accounted for, 
the zero value of roll amplitude occurs at a frequency other 
than resonant one of non-stabilized ship. It means that the ma- 
ximum amount of roll energy is absorbed by the tank from the 
ship at that other frequency. Hence in the case if sway is acco- 
unted for no phenomenon of double resonance can be directly 


observed. 


Influence of changing the tank parameters 
on course of roll transfer function 


It was assumed that in the preliminary design stage it is 
more convenient to make use of the full-form equations but 
without taking into account sway. Information this way obtai- 
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ned is more transparent and more useful for making relevant 
estimations and comparisons. 

In Fig.10 based on Fig.9 , only the courses for the dimen- 
sionless damping coefficients of liquid in tank, uoo = œ and 
Ugo = 0.8, are left. Fig.10 is further used asa a basis for making 
comparisons of tank stabilizing effectiveness depending on 
changes introduced to its location or other parameters. Fig.11 
shows in which way the roll reduction effect is influenced by 
a change in vertical location of the stabilizing tank. Next, in 
Fig.12 it is presented to which extent a change of tank meta- 
centric radius impairs ship roll amplitudes. From the presented 
diagrams it can be easily observed that their courses are almost 
in the same way unfavourably influenced by the decreasing of 


AA 
IZA 


—_——s 
———— 


0 
0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 
Fig.9. Diagrams of roll transfer functions without taking into account sway 


0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 
Fig. 10. Main transfer functions, without accounting for sway 
kg =0.97 ; € =0.013 ; 2 =20.00m ; rz =63.33m 


Ma 


Z 


0.40 0.60 0.80 1.00 1.20 1.40 
Fig. 11. Influence of the vertical coordinate of tank location, 
Zp=2.00m ; Zp = 20.00 m (Fig. 10) 


1.60 


tank elevation and the decreasing of tank metacentric radius. 
In both the cases the stabilizing effect is lower because the roll 
amplitudes of the stabilized ship are greater. By making use of 
the physical model (Fig.5, Fig.6, Fig.7) it can be observed that 
in both cases this is caused by changing the distance between 
the tank pendulum suspension point and the pivoting point of 
the whole system of pendulae of the physical model. 


a ppe 
ee i a R 
a 
EA A 
yy AN A 


0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 


Fig. 12. Influence of the tank metacentric radius - change of type of tank 
rz =46.73m ; rz = 63.33 m (Fig. 10) 


= Generally, it can be stated that the higher located the passi- 
ve stabilizing tank metacentre over the ship rolling axis, 
the greater the stabilizing effects. This is an important gu- 
ideline for designers. 


= Another important conclusion is : that type of selected pas- 
sive stabilizing tank is better whose metacentric radius is 
the greatest, at maintaining the remaining parameters, such 
as working liquid mass and its natural period of motion in 
tank, constant. 


= It is proposed to assume the ratio of the length of the tank 
metacentric radius and that of the tank mathematical pen- 
dulum as a relative measure of tank stabilizing quality, in- 
dependent on tank dimensions and mass of its working li- 
quid, and expressed as follows : 
Cini (25) 
For tanks of perpendicular sides the quantity C, is practi- 
cally constant, independent on their dimensions and the wor- 
king liquid level h, . For instance for the cubicoid tank and 
"Flume" tank it amounts to: Cz; = 0.82 and C72 = 0.67, respec- 
tively (the Cz2 value deals with the "Flume" tank having the 
dimensions of the middle throat equal to a half of tank length 
and a half of tank breadth, respectively). Cz values start to chan- 
ge distinctly for greater values of the liquid level in tank, h, , 
(i.e. for h, > 0.2b,) which is associated with the necessity of 
using more complex relationships in determining the natural 
periods of working liquid motions. Basing on the above given 
statements one can postulate the following thesis : 


Effectiveness of the passive stabilizing tanks of free liquid 
surface increases when the ratio of their metacentric radius 
and length of mathematical pendulum increases; 
the pendulum length is determined by the natural period of 
working liquid in tank. The ratio in question can be called 

"the stabilizing effectiveness of tank of a given type". 
The presented relationships open an interesting area for 
searching for optimum solutions of passive stabilizing tanks 
having free working liquid surface, 
and also for further research on ship roll stabilization. 


RECAPITULATION AND CONCLUSIONS 


O The presented mathematical model can be deemed well 
adjusted for designing the passive stabilizing tank to be in- 
stalled in ship. The elaborated calculation method can be 
used as a tool for aiding preliminary design stage of a ship 
fitted with passive stabilizing tank. 


O For correct design of stabilizing tanks the following three 
principles formulated on the basis of the presented model 
of ship with stabilizing tank, are important : 


** the higher located metacentre of passive stabilizing tank 
over ship rolling axis, the greater expected stabilizing 
effects 

** that type of passive stabilizing tank is better whose me- 
tacentric radius is the greatest, at maintaining the remai- 
ning parameters such as working liquid mass and natu- 
ral period of motions of liquid, constant 

** the higher the ratio of the tank metacentric radius and 

the mathematical pendulum length determined by na- 

tural period of motions of liquid in tank, the higher the 
stabilizing effectiveness of a selected type of passive 
stabilizing tank with free liquid surface. 


Moreover, it should be added that : 


> in spite of the simplifications resulting from the assu- 
med physical model, the achieved solutions provide suf- 
ficiently exact description of roll of a ship fitted with 
stabilizing tank. The wide range of the simplifications 
was introduced in order to reach "design usufulness" of 
the elaborated equations and their solutions 

> the way of accounting for sway in the obtained equa- 
tions arouses some reservations with regard to the in- 
troduced important simplifications. Nonetheless it was 
proved that in the case of real sea conditions and acco- 
unting for sway the so called double resonance did not 
directly appear for a ship with stabilizing tank. For this 
reason its stabilizing effect may appear lower than that 


expected. 
NOMENCLATURE 
aw - distance between rolling axis and waterplane, [m] 
D, D, - buoyance force of ship with and without tank, 
respectively, [kN] 
h, h, - initial metacentric height, not corrected regarding 


free-surface influence, for ship with and without 
tank, respectively, [m] 
ko, Uo and Ugo - ship roll damping coefficient, kg m2/s, 1/s, and [-] 
Ko» Ug and Ugo - damping coefficient of motions of liquid in tank, 
[kg m2/s, 1/s] and [-] 


ko = We/Wo - relative natural frequency of motions of liquid in 
tank, tuning factor, [-] 

lo - length of ship mathematical pendulum, [m] 

lip - length of tank mathematical pendulum, [m] 

M; - ship roll excitation moment, [Nm] 

Ms - ship mass, [t] 

m, - mass of liquid in tank, [t] 

mo - mass of ship mathematical pendulum, [t] 

Mo - mass of tank mathematical pendulum, [t] 

R - distance between metacentre of tank and 
metacentre of increment of ship immersed volume, 
corresponding to mass of liquid in tank, [m] 

R' - distance in physical model corresponding to R of 
ship, [m] 

Ro - distance between tank metacentre and ship rolling 
axis, [m] 

R'p - distance between tank pendulum suspension point 
and physical model pivoting axis, [m] 

ts - metacentric radius of ship without tank, [m] 
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- ship mass inertia radius, [m] 
- tank metacentric radius, [m] 


Tom - mass inertia radius of liquid in tank, [m] 

t - time, [s] 

T - draft, [m] 

Vi - horizontal component of velocity of mass of liquid 
in tank, due to ship sway, [m/s] 

Yn - ship sway, [m] 

ZG - height of mass centre of ship with tank, over 
reference plane, [m] 

ZGs - height of mass centre of ship without tank, over 


reference plane, [m] 

- height of mass centre of liquid in tank, over 
reference plane, [m] 

a - waterplane coefficient, [-] 

- amplitude and effective amplitude of wave slope 
angle, respectively, [rad, °] 

ô - ship block coefficent, [-] 


Kp, Ky - correction coefficients of wave slope angle, 
dependent on B/A and T/A, where : A - length of 
ship roll inducing wave 

E = m,m, - relative mass of liquid in tank, [-] 

Pe - metacentric radius of displacement layer of AT 
thickness, [m] 

To - natural oscillation period of ship pendulum and 
real ship, [s] 

To - natural oscillation period of tank pendulum and of 
liquid's motions in tank, [s] 

® - ship heeling angle, [rad, °] 

C0) - angular displacement of mass centre of liquid in 
tank, respective to ship, (mean slope angle of 
liquid surface in tank, respective to ship), [rad, °] 

= /M -relative excitation frequency, [-] 

(OJN - wave frequency, [1/s] 

Op - ship roll natural frequency, [1/s] 

Wg - natural frequency of liquid's motions in tank, [1/s] 
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m C onference 4 
Scientific Seminar 
E E 


of Regional Group 
of the Section on Exploitation Foundations 
On 24 June 2004 Gdynia Maritime University (GMU) 
hosted the 2nd — in-this-year scientific seminar of the 
Regional Group of the Section on Exploitation Founda- 


tions, Machine Building Committee, Polish Academy 
of Sciences. 


Four papers whose authors came from 
the scientific staff of GMU Mechanical Faculty, 
were presented during the seminar : 


* Results of service investigations of transverse sliding 
bearings lubricated with the use of non-Newtonian oils 
by A. Miszczak, D.Sc. 

* Problems associated with modeling NOx emission from 
ship two-stroke engine — by J. Kowalski, M.Sc. 

* Influence of lubricating oil contamination on wear of 
elements working in mixed friction conditions 
by A. Młynarczyk, M.Sc., Eng. 

* Selected problems associated with diagnostics of one- 
-stage refrigerating system — by T. Hajduk, M.Sc. 


The seminar was ended by presentation of the labora- 
tory facilities of GMU Mechanical Faculty to the seminar 
participants. 
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Fatigue “safe-life” criterion for metal 
elements under multiaxial static 
and dynamic random loads 


Janusz Kolenda 
Gdańsk University of Technology 
Naval University of Gdynia 


ABSTRACT 


Random stress with Cartesian components of known statistical moments of their mean 
values and power spectral densities of stochastic stress processes is considered. To acco- 
unt for the mean stress effect the generalised Soderberg criterion for ductile materials is 
employed. An equivalent stress with periodic (in the mean-square sense) components is 
defined by means of the equivalence conditions based on the average strain distortion 
energy. Also, is formulated the fatigue “safe-life” design criterion which covers the condi- 
tions of both static strength and fatigue safety and includes yield strengths and fatigue 


limits ne have simple physical interpretation, can be determined by uniaxial tests, are directly related 
to the applied loads, and can reflect material anisotropy. 


Key words : design criteria, multiaxial loading, random stress, mean stress effect 


INTRODUCTION 


The past decade has shown that fatigue is still a great chal- 
lenge for the engineering community. The reasons of it ,among 
other, may be the random or stochastic character of most loads 
that occur in nature. Therefore in this paper being direct conti- 
nuation of the author’s paper [1] devoted to fatigue safety of 
metal elements under deterministic loads, an attempt is made 
to formulate the fatigue “safe-life” design criterion for metal 
elements by using probabilistic approach. It is assumed that 
the considered stress is described by its Cartesian components 
as follows : 


9, (t)= Cj +s; (t) ? i= X, Y, Z, XY, YZ, ZX (1) 
where : 
cj - random mean values of known ot 
ments 
si(t) - zero mean stochastic processes of known power 
spectral densities, which are stationary (in the wide 


sense), stationary correlated with each other and 
statistically independent of the mean values c; . 


statistical mo- 


For the sake of brevity the stress components 
0,(t), Oy(t) and 6,(t) are dropped. 


In view of practical calculations, the actual stress compo- 
nents should be modelled with equivalent stress components 
of arelatively simpler form by means of an appropriate fatigue 
strength theory. For the multiaxial stress when axial forces, 
bending moments and torsional loads vary in time none of the 
fatigue strength theories is universally accepted and all fatigue 
criteria usually demonstrate large scatter [2]. On the other hand, 


in certain dynamic cases the conventional strength theories are 
considered to be satisfactory [3+5]. From some fatigue tests it 
has been concluded that also the criterion based on the average 
strain energy appears promising [6]. However, bearing in mind 
that for ductile metals the Huber-von Mises-Hencky strength 
theory is commonly accepted [7, 8], the average strain distor- 
tion energy is here, like in [1] , employed. 

The presented paper is aimed at the design criterion which 
would cover both the static strength conditions and fatigue sa- 
fety requirements under multiaxial stress. For this purpose use 
can be made of the generalised Soderberg criterion for non- 
-zero mean in-phase stress [1] : 


1/2 
py es |e 
i R; R,Ry 
1/2 (2) 
a xady f 
HDI <l , 1=X, y, Xy 
i F RF, 
where : 


Ri -yield strength relevant to the mean value of 
i-th stress component, ci 

F; -fatigue limit under fully reversed load relevant 
to the amplitude of i-th stress component, a; . 


As far as the computational effectiveness is concerned, spec- 
tral criteria can be very advantageous [9]. Therefore the ap- 
plied equivalence conditions are transformed into the frequen- 
cy domain. 
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EQUIVALENT STRESS 
UNDER MULTIAXIAL STATIC 
AND DYNAMIC LOADS 


The equation (2) suggests to model the stress 
components (1) by the equivalent stress components : 


ofe + $06) 


where : 


, 1=X,y,xy (3) 


- random mean value of i-th equivalent stress 
component 


s(t) - i-th zero mean stochastic process. 


It is assumed that the processes s(°4)(t) are periodically sta- 
tionary (in the mean-square sense [10]), stationary correlated 
with each other and statistically independent of the mean valu- 
es c’, which are sought in the form : 


sV(t)= a) sin(o,,t+@ ;)= (4) 
= ai exp(jogqt)+ aiz expl- jot) , 1=X, y, Xy 
where : 


ao, - random amplitude and phase angle 
of i-th equivalent stress component. 


(ea) i 
an = F exp jọ) > Ain = aiy (5) 
J 
(ai) = (ai) = (aian) = (anan) =0 (6) 
(iana) = (atay) =0 (7) 


Weg - equivalent circular frequency 


(e` - complex conjugate 
(*) 
In order to calculate the amplitudes of the equivalent stress 


components the following equivalence condition [1] is used at 
the beginning : 


- expected value 


1 T 1 T 
T Is cq (t)at = = fo (t)dt (8) 
where : T is the averaging time 
(eq) P (eq) ? 
oll + lo + 
ta (1) = ae, 9) 


3E | gg) +3 bF 


is the strain energy of distortion per unit volume 
in the equivalent stress state and 


TOEA s +04 =0,0, +302 ) (10) 


is that in the actual stress state. 


By applying (8) through (10), 
the following equivalence conditions can be written : 


+ [lem ofar= iOa. ix on 


22 POLISH MARITIME RESEARCH, No 3/2004 


ae : ji 
= [ol (oe (that =z Jox(o,(at (12) 
0 


Such equations of integral time averages are not conve- 
nient for evaluation of parameters of the equivalent stress com- 
ponents by using spectral data therefore they were replaced by 
the equations of ensemble averages : 


(bf) = (a? (t)) »1=X,y,XY (13) 
(206%) = (2.6020) 


To obtain the frequency domain formulation of (13) they 
were rewritten in terms of correlation functions as follows : 


(14) 


aj exp(-jouts}ralexplioggt N 
. Bo +a; expli Oeqt2 )+ aiz exp(- jJOeqts ) a 


= (fe; +sf t Jle: +s; t2 ))) 


Stationarity, in the wide sense, 
of the processes s;(t) implies that : 


(lc; +8; (tille; + silta J) = (c?)+ K,(t) a9 
where : 
Pa 
K ,(t)= (si (t )si(t2)) 


is the autocorrelation function of the process s;(t). In accor- 
dance with (3) through (6), and (16), one gets from (15) : 


(ef) (lf errlione) 


(17) 
+exp 2 j oT) = (o? ) +K,(t) 
Thus : 
(lee) =) a» 
“(iF )fexp(io.ge)+ exp jot = Kilt) (19) 


Fourier transformation of (19) gives : 


Hey) (0 ~ eq )+ 8(0 +Oeq ) =a (0) (20) 


where : 
5 - Dirac’s delta function 
S;(@) - power spectral density of the process s;(t). 


The mean-square value of the amplitude of i-th equivalent 
stress component can be estimated with the use of (20) by its 
integration over the whole frequency range, which yields : 


se) - fs, (o)do (21) 
and : 
(( (a) P) =28,(o)do (22) 


Similarly, stationary cross - correlation > 1/2 
of the processes 0,(t) and oy(t) requires that : Pe (ea) aleda ea) cos lo, — Ox ) (32) 
‘ d = e ee 
(o%(t,)o, (t2) =(c,c,)+K,y (t) 3) -| F FF, 
where : After expanding the functions f7" and f;' into Taylor se- 
K. (t)= (s; t,)s_(t ) ries around expected values of their arguments and retaining 
xy ) ü ( l ) vf 2) the linear terms, one gets : 
is the cross-correlation function 1/2 
of the processes sx(t) and sy(t). ((c a (c (ea), a 
By proceeding with (14) — in the same way (e) = > ae ee YO (33) 
as above — the following was obtained : i R; R,Ry 
eq), (eq FACCINE 
(e! oy n ay i: 1/2 
COPAAN 
‘ x a: eq), teq = 
{exp|i(o,-9, JJexp (io .4t)+ (24) (r,!\= a i 7 (al ay coslo, Px )) 
+ ex e i( — Jlex a jo 1} = i i F? EE 
PL J Py — Px JESP UJO eg : a (34) 
= (c,c, ) i Kyy (:) Hence the fatigue “safe-life” design criterion for anisotro- 
So: pic metal elements under multiaxial static and dynamic ran- 
CAE dom loads becomes as follows : 
(oe ey as) (25) i 1/2 
aplilo,- o.) pl (es) 
j + e 7 (35) 
ex O aT 
LE RONCI PA og =K_, (2) 12 
x y > X.Y 0 
4 + exp ile,- 9, 4 2 f 5 Silo) _ Sat) do <1 
: 2 
* exp z Jo at) (26) -o| i F; RF 
After Fourier transformation of (26) one obtains : CONCLUSIONS 
exp lj ( -ọ ) . O The fatigue “safe-life” design criterion which covers the 
Á . conditions of both static strength and fatigue safety of me- 
1 (a) fea) -3(o — Deg )+ tal elements under multiaxial static and dynamic random 
— {av a , = Syy lo) loads, was formulated. 
+exp |- J\Py— Px ) i i O The presented criterion includes material constants which : 
have simple physical interpretation 
‘ ne < pte phy Tp 
Alo eq ) (27) e can be determined by uniaxial tests 
. . e are directly related to the applied load 
where : Sx y (œ) is the cross power spectral density of the pro- e and can reflect material anisotropy. 
cesses O,(t) and Oy(t). Integration of (27) over the whole fre- 
quency range yields : NOMENCLATURE 
1 aj - amplitude of i-th component of the in-phase stress 


5 (ala cos lo,- Px ) = in (@)do (28) 


1.¢.° 


(ava cos lo,- Px ) =2 [Sey (o)do (29) 


FATIGUE “SAFE-LIFE” CRITERION 


With regard to (2) through (4), the criterion in question 
can be formulated as follows : 


(f,')+(fy')<1 (30) 
where : 
(eq) 2 cee (ea) 
i i xy (31) 


(i =X, y, Z, XY, YZ, ZX) 
aco) - random amplitude of i-th equivalent stress component 
aj}, aj2- quantities defined by (5) 
Ci - mean value of i-th component of the in-phase stress, 
random mean value of i-th component of the actual stress 
cka) - random mean value of i-th equivalent stress component 


E - Young modulus 
F; - fatigue limit under fully reversed load relevant to the 
stress amplitude a; 


fy, fs - partial safety factors acc. to [1] 

j - imaginary unit 

K; -autocorrelation function of the process s;(t) 

Kxy - cross-correlation function of the processes s,(t) and s,(t) 
R; -yield strength relevant to the mean stress value c; 

Si - time-variable part of i-th stress component 

slo) - time-variable part of i-th equivalent stress component 
S; - power spectral density of the process s;(t) 

Sxy ~- cross power spectral density of the processes s,(t) and sy(t) 
t - time 


POLISH MARITIME RESEARCH, No 3/2004 23 


MARINE ENGINEERING 


T - averaging time 
ô - Dirac’s delta function 
v - Poisson’s ratio 
©; -i-th stress component 


os) - i-th equivalent stress component 
- time interval 


T 
@;  - random phase angle of i-th equivalent stress component 


, eq - strain energy of distortion per unit volume in the actual 
and equivalent stress states 

oO - circular frequency 

@eq - equivalent circular frequency 

(e) - expected value 


Q - complex conjugate 
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Miscellanea 
TOP KORAB Ministerial award 
activity in 2003 Some years ago a team of the Faculty of Ocean Engine- 
ering and Ship Technology, Gdansk University of Technolo- 
In 2003 The Polish Society of Naval Architects and 8Y. has commenced working on design and construction of 
Marine Engineers, TOP KORAB, — in accordance with its 4 Prototype fishing cutter for Baltic Sea, which has been ho- 
custom of arranging topical meetings — carried out, in ped to be useful in renewal of the obsolete fishing vessels 
Gdansk, 9 meetings devoted to the following topics : flying Polish flag. The undertaken task was sucessfully com- 
7 f — : pleted and the built cutter has already operated for a year, 
s Achievements of POL-LEVANT shipping company in gathering flattering opinions from the side of its users. 
, the light of situation of Polish shipping © Polish state authorities received with recognition that 
= Development stategy of Pomeranian region important achievement seeing in it a nucleus of a new gene- 
= Certification of management systems ration of ecolgical family-operated fishing vessels. Giving 
K Dimension allowance system applied in shipbuilding voice to it, the authorities rewarded the team consisted of the 
ae i it i 7 following persons : 
“* Current problems of Central Maritime Museum in . 
Gdansk - presentation of permanent exhibitions f - EE, 4 R Eng. 
“+ Achievements and development prospects of Gdynia AJ ite DE cae ng. 
Naval Shipyard A C D = D ri 8 E 
œ Prospects and future of Polish Shibuilding Industry Fo- A Gk R N Daa 
ium fi T. Blekiewicz, Eng. 
% A proposal of topics of club meetings to be organized in A E o 
2004, and of forms of future activities of TOP KORAB A M S ar 4 m ng. 
“+ Presentation of the current state of TOP KORAB chro- i ai tac owia a nr 
nicle. with the Ministerial prize which was solemnly 
handed over on 13 May 2004. 
A full-day coach excursion to Malbork to visit the his- i 
torical Teutonian Knights' castle, was also carried out. 
Whereas in Szczecin, apart from 4 organizational 
meetings, 2 topical meetings were held : 
* Current problems of NOWA Szczecin Shipyard Co 
* Problems associated with implementing a new IMO code: 
ISPS Code 2003. 
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Physical aspects of application and usefulness 
of semi-Markovian processes for modelling 
the processes occurring in operational 
phase of technical objects 


Jerzy Girtler 
Gdańsk University of Technology 


ABSTRACT 


In the paper usefulness of semi-Markovian processes for modelling real processes which 
occur in the operational phase of various technical objects, is considered. The usefulness 
of this theory was proved by presenting specific features of the semi-Markovian processes 
and physical aspects of their use as models of the processes occurring in the operational 
phase of the objects in question. The specific features were described with taking into 
account the process of changing the technical states of such objects. The physical aspects 
of models in the form of semi-Markovian processes were justified by postulating a hypo- 


thesis by which wear process of tribological units of crucial sub-assemblies of machines, can be explained. 
Also, the consequences logically resulting from the hypothesis and necessary to its verification, are discus- 
sed. A method for its verification is also attached. It was proved possible to model the real processes in 
question by means of the continuous semi-Markovian processes of finite set of states. 


Key words : technical object, semi-Markovian process, tribological unit, wear 


INTRODUCTION 


The theory of semi-Markovian processes gains more and 
more widespread applications in engineering sciences and in 
operation of technical objects (e.g. diesel engines, gas turbi- 
nes, screw propellers, pumps, compressors, coolers, filters, ship 
propulsion systems). It makes it possible to elaborate models 
of various real processes, including processes of changes of 
technical and operational states as well as operational proces- 
ses of any technical objects. 

The models are elaborated in the form of special stochastic 
processes, 1.e. semi - Markovian ones, and recently also deci- 
sion-related (controlled) semi-Markovian processes whose re- 
alizations depend on decisions made in the instants of chan- 
ging their states [4, 6, 8, 11]. However it is not easy to apply 
the processes because of their specific features. Not taking them 
into account may result in building such models of real proces- 
ses in the form of semi-Markovian ones whose investigations 
cannot provide any new information about a modelled pro- 
cess, regarding e.g. durability and reliability of technical ob- 
jects, their load spectra etc. Therefore it is worthwhile to indi- 
cate the physical aspects of application of semi-Markovian pro- 
cesses as models of real processes related to technical objects, 
at least with taking as an example the process of changing tech- 
nical states during operation of such objects. 


SPECIFIC FEATURES 
OF SEMI-MARKOVIAN PROCESSES 


The semi-Markovian processes are stochastical ones of 
peculiar features. In the literature there are different definitions 
of semi-Markovian process, which have different ranges of 
generality and exactness. For purposes of the modelling of ope- 
ration of technical objects the semi-Markovian process (fami- 
ly of random variables) 

{Y(t): te T} at T=[0, +] 
can be defined by means of the so - called uniform Markovian 
renewal process. Such definition proposed by F. Grabski [9] is 
close to those given by other authors [13, 19, 20, 1]. 

From the definition it results that it is stochastic one of 
a discrete set of states and its realizations are right-hand conti- 
nuous functions constant within intervals (of uniform values 
within operational time intervals which are random variables). 
Such process is defined only when its initial distribution P; = 
= P{Y(0) = si} as well as the functional matrix Q(t) = [Qi] is 
known; the matrix elements are the probabilities of transition 
from the state si to the state sj, within the time not greater than 
t #j;1,j = 1, 2,..., K), being the non-decreasing functions 
Qij(t) of variable t, namely : 


Og = PUY) ESj > ta Tn < t|Y(t,) =s;} (0) 
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A semi-Markovian model of an arbitrary real process can 
be formed only when states of the process can be defined in 
such a way that duration time of a state appearing in the instant 
Tn as Well as a state possible to be reached in the instant T,+1 do 
not depend stochastically on the preceding states and their du- 
ration time intervals. 

For elaborating the semi-Markovian model {W(t) : te T} 
of a given real process it is necessary to apply the theory of 
semi-Markovian processes. It makes it possible to determine 
probabilistic characteristics of an arbitrary random process 
(if only such model is formed), which may be of a practical 
importance. Such process may be a real process of changing 
technical states of any technical object. An example model of 
such process can be presented in the form of a graph of chan- 
ging the states of the object (Fig.1). The model is one of the 
crucial components of the operational model of every technical 
object. One of the possible realizations of the process of chan- 
ging the technical states, {W(t) : t€ T}, is presented in Fig. 2. 


T33 > P23 


Fig. 1. A graph of changing the states of an arbitrary technical object : 
s] - full serviceability state, s2 - partial serviceability state, 
s3 - unserviceability state, pij - probability of passing the process from the 
state si to the state s;, Tij - duration time of the state s; provided 
the process passed to the state Sj; ij=1, 2,3 


W(t) 


t 
To Ty T3 T3 TA T5 


Fig. 2. An example realization of the process {W(t) : t € T} of an arbitrary 
technical object : {W(t): t € T} - process of changing the technical states, 
t - time of operation; s] - state of full serviceability, s3 - state of partial 
serviceability, s3 - state of unserviceability 

In reality the process of changing the technical states of 
any technical object is stochastic, continuous in time and thro- 
ughout the states. It means that the realizable kinds of techni- 
cal state form a set of technical states of a given object, which 
is infinite one. 

Identification of all technical states of many technical ob- 
jects (e.g. combustion engines, boilers, pumps, compressors, 
ship prpulsion systems) is neither possible nor purposeful, both 
due to technical and economical reasons. Hence it is necessary 
to split the set of states of technical objects into a few number 
of classes (sub-sets) of technical states. 

Assuming, as a splitting criterion, usefulness ofa given tech- 
nical object for realization of tasks (serviceability) one can di- 
stinguish the following classes (sub-sets) of technical states, 
shortly called "states" [4, 6, 8] : 
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* the state of full serviceability, sı , which makes it possible 
to use the object in any conditions and in any range of lo- 
ads, to which it was adjusted in the phases of its designing 
and manufacturing 


* the state of partial serviceability, s, which makes it possi- 
ble to use the object in limited conditions and in a range of 
loads lower than those to which it was adjusted in the pha- 
ses of its designing and manufacturing 


* the state of unserviceability, s3, which does not make it 
possible to use the object in accordance with the purpose 
for which it was intended (e.g due to its failure, carrying 
out maintenance operations on its subassemblies etc.). 


Therefore it is a three-state process 
of continuous realizations (continuous with time). 


The technical objects which are in the full serviceability 
state (s1) can be used in any instant and conditions to which it 
was adjusted in the designing and manufacturing phases, and 
under various loads. And, the technical objects which are in 
the partial serviceability state (s2) can be used or maintained 
depending on a decision-making situation (control strategy), 
whereas the unserviceable technical objects (in the state s3) 
due to their failure are always maintained provided it would be 
cost-effective. However the technical objects unserviceable 
because of realization of preventive maintenance, are unservi- 
ceable only during the time of realization of these operations 
which require the object's structure to be trespassed by disas- 
sembling some of their devices. 


The particular states s; E€ S (i= 1, 2, 3) can be identified by 
means of an appropriate diagnostic system (SD) whose useful- 
ness depends on quality of an applied diagnosing system (SDG) 
as well as on its capability of identifying the states of the tech- 
nical object considered as a diagnosed system (SDN). 


It can be assumed that if whichever of the states s or s3 
does not occur then the object in question is in the state s4. 


The considered process of changing the technical states of 
technical object is, in mathematical terms, a function which 
maps the set of the instants T into the set of technical states, S. 
Such process can be modelled by means of the stochastic pro- 
cesses of discrete set of states and continuous duration time of 
distinguished technical states of the object. 


Therefore the set of technical states : 


S = {S1, 82, 83} 


can be considered as the set of values 
of the stochastic process 


{W(t): te T} 
whose realizations are constant within time intervals 
and right-hand-side continuous, Fig.2. 


In the case of the process of operation of technical objects 
the following characteristics may be of a practical importance : 


> the one-dimensional distribution of the process (instanta- 
neous distribution), whose elements are the functions P,(t) 
representing the probability of the event that in the instant t 
the process will enter the state sę 


> the limiting distribution of the process P ;= lim P{Y(t) = s j } 
tec 


> the conditional probabilities, i.e. those of transition of the 
process from the state s; to the state sj , 
Pij(t) = PL Y(t) = sj / Y(O) = si} 
(transition probabilities) 


> the distribution of the time of the first transition of the pro- 
cess from the state s; to the sub-set of the states A (®;,(t)), 
and if this sub-set contains only one element — to the state 
sj, i.e. the distribution ®;;(t) 


> the distribution of return time of the process to the state s;, 
i.e. the distribution ®;;(t) 


> the asymptotic distribution of the renewal process {V;,(t) : 
t 2 0}, generated by return time intervals of semi-Marko- 
vian process (to the state sj available from the state s;), which 
at the instant t, obtains a value equal to number of "co- 
ming-in" events of that process to the state sj 


> an approximate distribution of the total time of maintaining 
the process Y(t) in the state s; provided the state s; is that 
initial one 


> the expected value E(T;) of the duration time T; of the state 
s; of the process irrespective of the state to which transition 
occurs at the instant Tp+1 


> the variance D*(T;) of the duration time T; of the state s; 


> the expected value E(T;j) of the duration time Tj of the state 
s; of the process provided the state sj is the next one 


> the expected value E(O;;) of the random variable ©;j which 
represents the return time of the process to the state sj 


> the expected value E{V;;(t)} of the random variable V;,(t) 
which represents the number of "coming-in" events of that 
process to the state sj within the time interval [0 , t] 


> the variance D*{Vj(t)} of the random variable V(t) 


> the average number of "coming-in" events, À;j(t), of the pro- 
cess to the state s;, related to a unit of time, provided the 
state s; of the process is that initial one (i.e. the intensity of 
the " coming-in" events of the process to the state sj provi- 
ded Y(0) = s;) 


> the limiting intensity of the "coming-in" events of the pro- 
cess to the state s;, i.e. the intensity hi= limA ;,(t) ; 
t—e0 


To obtain numerical values of the above mentioned 
characteristics is possible if two following 
conditions will be satisfied : 


e ifappropriate statistical data whose values would represent 
estimation of the transition probability p;j, of the expected 
value E(T;), etc, will be collected 


e ifasemi-Markovian model of operation process of techni- 
cal objects having a small number of its states and mathe- 
matically simple functional matrix Q(t) will be elaborated. 


The second condition is important in the case of calculating 
the instantaneous distribution of states of the process P;(t). The 
distribution can be calculated if the initial distribution of the 
process and its functions Pjj(t) are known. The calculating of 
the probabilities P;j(t) consists in solving the set of Volterra se- 
cond-kind equations in which the functions Q;j(t) being elements 
of the process functional matrix Q(t), are known quantities [9]. 

In the case when number of process states is small and the 
functional matrix of the process — rather simple, that set can be 
solved by using the Laplace transform method [9, 10, 20]. How- 
ever when number of process states is large or when its func- 
tional matrix (core of the process) is very complex, only an 
approximate solution of the set of the equations is available. 
Such (numerical) solution does not make it possible to deter- 
mine values of probabilities of occurrence of process particu- 
lar states if t is of a large value (theoretically if t 4 ~). The 
numerical solution does not provide any answer to the question 


very important for operational practice, namely : in which way 
do the probabilities of semi-Markovian process states change 
if t is large? 

From the semi-Markovian process theory it results that the 
probabilities, in the case of the ergodic semi-Markovian pro- 
cesses, tend along with time to strictly determined constant 
numbers. They are called the limiting probabilities of states 
and their sequence forms the limiting distribution of the pro- 
cess. The distribution makes it possible to define the availability 
factor of technical object as well as the income or cost per unit 
time of its operation. The quantities serve as criterion func- 
tions in solving problems of operation process optimization of 
technical objects. Such distribution can be calculated much 
easier than the instantaneous one. 

Similar difficulties are associated with solving the set of 
integral equations which make it possible to calculate the di- 
stribution of the random variable Oja which determines the time 
of the first transition of the process from the state sj; to the sub- 
-set of states, A, or (if the state A contains only one element) — 
— to the random variable ©; which determines the time of the 
first transition of the process from the state s; to the state sj. 


PHYSICAL ASPECTS OF APPLICATION 
OF SEMI-MARKOVIAN PROCESSES 
FOR MODELING THE PROCESSES 
OCCURRING IN THE PHASE 
OF OPERATION 


From the presented considerations it results 
that the semi-Markovian models are characteristic 
of the following features [7, 9, 10, 19, 20] : 


¢ Firstly, is satisfied the Markov condition that future evolu- 
tion of the investigated object (e.g. the process of changing 
the technical states in the object's operational phase), for 
which the semi-Markovian model has been elaborated, sho- 
uld depend only on its state at a given instant but not on the 
functioning of the object in the future, i.e. that the future of 
the object would not depend on its history but only on its 
present. 


@ Secondly, the random variables T; (determining duration 
time of the state s; regardless which state will occur after it) 
as well as Tij (determining duration time of the state s; pro- 
vided the next state of the process will be the state s;) have 
their distributions different from exponential ones. 


Therefore in the modeling aimed at elaboration of a semi- 
-Markovian model of the process of changing the object's tech- 
nical states, an analysis of changing the states of real process, 
i.e. those occurring in the operational phase of the object in 
question, should be taken into account. 

In the case of every technical object the process of chan- 
ging its technical states is that where the duration time inter- 
vals of each its state are random variables. Particular realiza- 
tions of the random variables depend on many factors, a.o. on 
the technical object's wear. 

In the case of such technical objects as e.g. diesel engines, 
compressors or pumps, it was observed that wear of their sli- 
ding tribological units is weakly correlated with time [3, 8, 16, 
17, 23, 24]. The observation is important because serviceabili- 
ty of such machines depends mainly on the technical state (1.e. 
on wear) of their tribological units. This made it possible to 
predict technical states of such machines with taking into ac- 
count solely their present state and neglecting those earlier oc- 
curred. An explanation of the fact would make it possible to 
elaborate (by applying the theory of semi-Markovian proces- 
ses) more adequate mathematical probabilistic models for pre- 
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dicting the technical states of particular machines. To this end 
the following hypothesis (H) can be offerred : 


A state of an arbitrary sliding tribological unit as well as 
its duration time essentially depend on the state preceding it 
and neither on the earlier occurred ones and nor their dura- 
tion time intervals because its load and both rate and incre- 
ments of wear, implicated by it, are the processes of asymptoti- 
cally independent values. 


The last statement of the hypothesis (because its load...) 
results from two obvious facts : 


© there is a strict relationship between loading on sliding tri- 
bological units and their wear [15, 17, 23] 


= there is a lack of monotonic changes of loading on tribolo- 
gical units of machines within longer time of their opera- 
tion, hence their loading can be assumed stationary [3, 17, 
21,22,23]: 


The load stationarity (in a broader sense) means in every 
case that all multi-dimensional probability density functions 
depend only on mutual distances of the instants T1, T2,..., Tn , 
and not on their values [5]. Therefore the one-dimensional pro- 
bability density function of load values does not depend on the 
instant related to a given value, and the two-dimensional pro- 
bability density function depends only on difference of the in- 
stants in which observed loading values occurr. 

And, in a narrower sense, the fully stationary loading is 
understood as that whose all possible statistical moments of 
higher orders as well as the total moments of loading conside- 
red as a process, are not time - dependent. In the case of the 
fully stationary process (in a narrower sense) its characteristic 
quantities are as follows : 


the expected value m(t) =m = const 
variance D(t) = o^ = const 
autocorrelation A(T1 , T2) = A*(T2 - T1) = A*(r) 
and autovariance K(T1 , T2) = K*(t2 - T1) = K*(r). 


However the stationary process in a broader sense 
is characteristic of : 


m(t) =m = const as well as A(T, T2) = A*(T2 - T1) = A*(r). 


In practice the loading stationarity in a broader sense is 
important. And, in this case to investigate the loading on tribo- 
logical units in order to reveal the enumerated properties, is 
not necessary, as it is known from investigations of different 
machines, which have been performed so far, that the loading 
on their tribological units continously changes in such a way 
that its particular values measured after very short time inter- 
vals are strongly correlated to each other. However when the 
time interval between measurements of loads increases, the 
correlation between the loads decreases. Therefore the loading 
values measured in the time intervals (or instants) very distant 
apart can be considered as independent ones. This feature is 
called the asymptotic independence of a loading value measu- 
red in the instant, e.g. T+), from that measured in the instant Ti, 
when the range At = Ti+; - Ti is large enough. 

The so understood asymptotic independence between lo- 
ading values either measured or calculated in the instants 7; 
and Ti+1, is manifested by that their mutual dependence decre- 
ases along with increasing the range At. Moreover, from work 
principles of particular machines it is also known that their 
loading considered within a longer time of their correct opera- 
tion, does not (and cannot) show any monotonically increasing 
or decreasing changes. Therefore one can assume that the ma- 
ximum loading values appear in the specified instants acciden- 
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tally, always with some probability only. This lack of monoto- 
nicity of the loading is called its stationarity. 

In order to verify the presented hypothesis (H) it is neces- 
sary to predict the consequences whose occurrence can be con- 
firmed empirically if the hypothesis is true. The consequen- 
ces (K) which can be derived from the hypothesis (with ta- 
king into account the mentioned features of loading on ship 
power plant machines and their sliding tribological units) are 
the following [3] : 


ve K; - irregular course of realization of wear process of par- 
ticular sliding tribological units 


ve K, -interweaving realzations of wear processes of sliding 
tribological units 


v« K; - such course of autocorrelation function for a given 
sliding tribological unit that at first the function fast decre- 
ases along with the range 0 = hAt (h = 1, 2,..., n) incre- 
asing, and next it oscillates about zero, at relatively small 
amplitude smaller and smaller along with At increasing 


ve K; - almost normal distribution of wear increments of sli- 
ding tribological units for a sufficiently long time interval 
(At) of their correct operation 


ve Ks - linear relationship of variance of wear process of sli- 
ding tribological units and their operation time values. 


The above described consequences 
are graphically illustrated in Fig. 3 + 6. 


The presented consequences can be justified as follows : 


If the features of the loading on machines and thus their 
tribological units are such as those above described, then the 
course of realization of wear of the units will have to be irregu- 
lar. This forms the basis to assume that the wear increments 
recorded in the time intervals much distant from each other, 


z(t) 


Zlim 


t 
Fig. 3. Example realizations of wear processes of sliding tribological units : 
Z - WEAF, Zļim - limiting wear value, t - object's operation time 


r(0) 
K3 


1,0 


2 


0 
Fig. 4. Example course of the autocorrelation function r(0) where : 
0- the range between time intervals (within which wear was investigated) 


are asymptotically independent, and along with increasing the 
time (time range 0, where : 0 = hAt, h = 1, 2,..., n) between the 
intervals the relationship between the increments in question 
will be weakening. Hence the wear processes of such units can 
be considered as those of asymptotically independent incre- 
ments [3]. 


Au 


EM) 
Fig. 5. Example form of the density function f (Au) of the asymptotically 
independent wear increments Au of sliding tribological units, 
for sufficiently long time interval (At) of their correct operation 


D{Z(t)} Ks 


t 
Fig. 6. Example relationship of the variance DZ) )} of the wear process Z(t) 
of sliding tribological units of asymptotically independent increments, 
and their operation time t 


One of the most important characteristics of every stocha- 
stic process (including wear processes) is the autocorrelation 
function r(0) = f(0), where : r - autocorrelation coefficient, 
O - time range. If r(0) decreases then wear process of a given 
tribological unit can be considered as that of asymptotically 
independent increments. Hence for tribological units in question 
an almost normal distribution of wear increments should be 
expected in a sufficiently long time interval. Moreover, if wear 
process is of asymptotically independent increments its variance 
D {Z(t)} increases linearly along with correct operation time, 
in accordance with the relationship (2) : 


D?{Z()} = At +B 


where : 
A and B - process constants 


(2) 


The above mentioned consequences K;(i = 1, 2,..., 5) reve- 
al the probabilistic principle of wear of sliding tribological units. 
They are not contradictory to each other, and their logical ve- 
racity is doubtless. Hence it is possible to consider all of them 
as one consequence K and to use it for empirical proving if the 
presented hypothesis (H) is true or false. Such verification con- 
sists in experimental testing of wear of sliding tribological units 
and checking if the consequence K is true, which is equivalent 
to checking if the consequences K; (the facts) appear or do not 
appear. Such verification of the hypothesis H makes it neces- 
sary to accept the veracity of the following syntactic implica- 
tion [8, 12, 18]: 


H>K (3) 


Then the non-deductive (inductive) reasoning may be used 
in accordance with the following scheme [8 , 18] : 


(K,H=>K)tH (4) 


where : 
K= {K;,1=1, 2,..., 5} 


Its logical interpretation is as follows : if the empirical te- 
sting of the consequence K confirmed its veracity, then if the 
implication (3)is true, the hypothesis H is also true and accep- 
table. The reasoning in accordance with (4), called reductive 
one, does not lead to any firm conclusions but only to probable 
ones [12 , 18]. 

Semi-Markovian model of an arbitrary operational process 
can be applied in the case of machines even by using the dia- 
gnostic systems (SDG) where the reductive reasoning can be 
applied. Hence when reasoning on a diagnosis (conclusion) 
concerning the state of a given technical object being a diagno- 
sed system (SDN) the statement K (stating that this — and not 
another — vector of values of diagnostic parameters is obser- 
ved) is deemed to be a fully firm premise. 

However the statement S (stating that there is this — and not 
another — state of SDN) is a conclusion formulated on the basis 
of the statement K during the non-deductive reasoning process 
which proceeds in compliance with the following scheme : 


(K,S=>K)tS 


where : 


(5) 


K - fully firm premise 
S - conclusion formulated on 
the basis of the statement K. 


From such reasoning the following hypothesis results : 
The considered SDN is in the state S because the vector 
of values of diagnostic parameters, K, is observed. 


It can be also formulated in another, equivalent way : 
The vector of values of diagnostic parameters, K, 
is observed because the SDN in question is in the state S. 


Such reasoning does not make it possible to formulate firm 
conclusions but only probable ones. Therefore it is not possi- 
ble to exactly determine a technical state of a diagnosed sys- 
tem (SDN) and thus to control its operational process in such 
a way that a future state would depend on many previous states. 


FINAL REMARKS AND CONCLUSIONS 


% The semi-Markovian processes are useful models for inve- 
stigating the real processes which occur in the operational 
phase of technical objects. Hence elaboration of a semi- 
-Markovian model of a given process occurring in opera- 
tional phase of an arbitrary technical object, makes it pos- 
sible to easily determine probabilistic characteristics of the 
process in question. 


** In practice the semi-Markovian processes are more useful 
than Markovian ones. It means that the semi-Markovian 
processes of continuous time parameter and finite set of 
states are characteristic of that the time intervals of staying 
the processes in particular states are random variables of 
arbitrary distributions concentrated in the set R+ = [0 , œ]. 
This differs them from the Markovian processes whose in- 
tervals are random variables of exponential distributions. 


* A semi-Markovian model of an arbitrary process occurring 
in the phase of operation of a technical object is that of 
a finite set of states and continuous time. 


s% An additional benefit from application of semi-Markovian 
processes (like also of Markovian ones) is the possibility of 
using professional computer tools for solving different sets 
of state equations for such models of real processes. 
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Semi-Markovian model of an arbitrary operational process 
can be applied in the case of machines even by using the 
diagnostic systems (SDG) where the reductive reasoning 
can be applied. However such reasoning does not make it 
possible to formulate firm conclusions but only probable 
ones. 
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Plenary session of Machine Building Committee 
held in coast region 


On 31 May and 1 June 2004 members of the Machine 
Building Committee, Polish Academy of Sciences held, this 
turn in Gdansk, their plenary session commencing the next 
tenure of activity of the Committee. It was devoted, apart 
from organizational matters, to new challenges arising from 
the entrance of Poland into European Union. 


On this occasion took place an open scientific meeting on : 


Development prospects of machine building 
and operation after entrance 
of Poland into European Union 


8 papers were prepared to 
be presented during the meeting : 


> Poland in the European Space of Science and Education 
by Prof. W. Sadowski (Gdansk University of Technology) 

> Polish Space of Science — prospects of development 
by Prof. J. Kicinski (Insitute of Fluid Flow Machinery, 
Polish Academy of Sciences, Gdansk) 

> Application trends in EU programs in the area of machi- 
ne building and operation — by Prof. A. Mazurkiewicz 


C (Institute of Technical Operation Processes) 
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> Contemporary research and education problems of ocean 
engineering — towards European Space of Scientific Rese- 
arch — by Prof. J. Szantyr (Faculty of Ocean Engineering 
and Ship Technology, Gdansk University of Technology) 

> Proecological activity of Mechanical Faculty concerning 
machine building and operation, in the frame of EU 
by Prof. W. Przybylski (Mechanical Faculty, Gdansk Uni- 
versity of Technology) 

> Attempts to problems of operation of ship engines in the 
aspect of cooperation in the frames of EU and NATO 
by Prof. R. Cwilewicz (Gdynia Maritime University), 
Prof. L. Piaseczny (Polish Naval University, Gdynia) 

> Integration of doctorate studies — by Prof. B. Żółtowski 
(Technical Agricultural Academy, Bydgoszcz) 

> MARIE CURIE individual grants — basic information 
by Prof. W. Zwierzycki (Poznan University of Techno- 


logy) 


The two-day meeting was co-organized by the Institute 
of Fluid Flow Machinery, Polish Academy of Sciences, Me- 
chanical Faculty and Faculty of Ocean Engineering and Ship 
Technology, Gdansk University of Technology. a 
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Miscellanea 


Reliability of Large Systems 


An interesting book written by : 


Krzysztof Kolowrocki 
Gdynia Maritime University, 


has been recently published by Elsevier (April 2004). 


The Book Website is : http://www.elsevier.com/wps/find/ 
bookdescription.cws_home/702918/description 


Reliability of Large Systems 


Krzysztof Kolowrocki 


The book consists of eight chapters 
whose content is shortly characterized below. 


I. Basic notions 


Basic notions and agreements which are necessary to fur- 
ther considerations, are introduced. The asymptotic appro- 
ach to the system reliability investigation and the system 
limit reliability function is defined. 


2. Two-state systems 


Two-state homogeneous and non-homogeneous series, 
parallel, "m out ofn", series-parallel and parallel-series sys- 
tems are defined. Their exact reliability functions are deter- 
mined. 


3. Multi-state systems 


Basic notions of the system multi-state reliability analy- 
sis are introduced. The multi-state homogeneous and non- 
-homogeneous series, parallel, "m out of n", series-parallel 
and parallel-series systems with degrading components are 
defined and their exact reliability functions are determined. 
The multi-state limit reliability function of the system, its 


risk function and other multi-state system reliability charac- 
teristics are introduced and determined. 


4. Reliability of large two-state systems 


Auxiliary theorems on limit reliability functions of large 
two-state systems, which are necessary for their approxima- 
te reliability evaluation, are formulated. The classes of limit 
reliability functions for homogeneous and non-homogene- 
ous series, parallel, series-parallel and parallel-series systems 
and for a homogeneous "m out of n" system are established. 

Applications of the asymptotic approach to reliability eva- 
luations of model systems are presented. Six corollaries are 
formulated and proved on the basis of the auxiliary theorems 
and applied to finding limit reliability functions of the con- 
sidered systems and approximate evaluations of their relia- 
bility functions, lifetime mean values and lifetime standard 
deviations. 

The reliability evaluation is done for the following sys- 
tems: the model non-homogeneous series system, the ho- 
mogeneous parallel system of an energy cable, the "16 out 
of 35" lighting system, the homogeneous regular series-pa- 
rallel gas distribution system, the non-homogeneous regular 
series-parallel water supply system and the model homoge- 
neous regular parallel-series system. 

The accuracy of the performed evaluations is illustrated 
in tables and figures. The reliability data of components are 
assumed either to be arbitrary or to come from experts. The- 
se reliability evaluations of the considered systems’ charac- 
teristics are an illustration of the possibility of applying the 
asymptotic approach to system reliability analysis of large 
real technical systems. 


5. Reliability of large multi-state systems 


Auxiliary theorems on limit reliability functions of mul- 
ti-state systems, which are necessary for their approximate 
reliability evaluation, are formulated and proved. The clas- 
ses of limit reliability functions for homogeneous and non- 
-homogeneous series, parallel, series-parallel and parallel- 
-series multi-state systems and for a homogeneous multi- 
state "m out ofn" system are established. Practical applica- 
tions of the multi-state asymptotic approach to reliability eva- 
luation of real technical systems are presented. 

On the basis of auxiliary theorems some corollaries are 
formulated and proved and then applied to approximate re- 
liability and risk characteristics determination of real techni- 
cal multi-state systems having series, parallel, "m out of n", 
series-parallel and parallel-series reliability structures. Eva- 
luations are given of multi-state reliability functions, mean 
sojourn times in the state subsets and their standard devia- 
tions, mean lifetimes in the states, risk functions, and exce- 
eding moments of a permitted risk level for selected real 
systems. 

The homogeneous series piping transportation system, 
the model homogeneous series telecommunication network, 
the homogeneous series bus transportation system, the non- 
-homogeneous series piping transportation system, the ho- 


POLISH MARITIME RESEARCH, No 3/2004 


31 


mogeneous parallel system of an electrical cable, the non- 
-homogeneous parallel rope system, the "10 out of 36" homo- 
geneous steel rope system, the model homogeneous series-pa- 
rallel system, the homogeneous and non-homogeneous series- 
-parallel pipeline systems and the homogeneous parallel-se- 
ries electrical energy distribution system are analysed and the- 
ir reliability characteristics are evaluated. 

Necessary data on system components reliability and sys- 
tem operation processes come from experts, from trade stan- 
dards and from certificates issued by the system component 
producers. Component reliability and system operation pro- 
cesses data are, out of necessity, approximate and they concern 
the components' mean lifetimes in the reliability state subsets 
and the hypothetical distributions of these lifetimes. The accu- 
racy of the asymptotic approach to the reliability evaluation of 
the considered systems is illustrated in tables and figures. 


6. Reliability evaluation of port 
and shipyard transportation systems 


The multi-state asymptotic approach is applied to the relia- 
bility and risk characteristics evaluation of selected large trans- 
portation systems used in ports and shipyards. Reliability ana- 
lysis of multi-state series, series-parallel and parallel-series 
transportation systems is based on corollaries formulated and 
proved in this and the preceding chapters. The corollaries are 
applied to evaluate reliability characteristics of three transpor- 
tation systems used at the Port of Gdynia and one operating at 
the Naval Shipyard of Gdynia. 

The port grain transportation system built of three-state non- 
-homogeneous series-parallel subsystems, the port oil piping 
transportation system composed of three-state non-homoge- 
neous series-parallel subsystems, the port bulk transportation 
system built of four-state non-homogeneous series-parallel and 
series subsystems and the shipyard rope transportation system 
that is a four-state homogeneous parallel-series system, are 
considered. 

Multi-state reliability functions, mean values of sojourn ti- 
mes in the state subsets and their standard deviations, mean 
values of lifetimes in the particular states, risk functions, and 
exceeding moments of a permitted risk level are determined 
for these systems. 

The accuracy of the asymptotic approach to the reliability 
evaluation of these systems is also illustrated. System compo- 
nents reliability data and system operation processes data come 
from operators of these systems, component technical certifi- 
cates and binding standards. Reliability data are, out of neces- 
sity, approximate and they concern only the mean values of the 
system components' sojourn times in the state subsets and hy- 
pothetical distributions of these lifetimes. 


7. Reliability of large multi-state 
exponential systems 


Limit reliability functions of multi-state series, parallel, 
"m out of n", series-parallel and parallel-series systems com- 
posed of components having exponential reliability functions, 
are established. Next, the results are presented in the form 
of tables containing exact algorithms of the procedure while 
evaluating reliability characteristics of these systems' reliabili- 
ty in order to deliver a simple and convenient tool for everyday 
practice of the reliability practitioners. The tables are compo- 
sed of three parts, containing reliability data of the evaluated 
system, necessary calculations, and results of the system relia- 
bility evaluation. The way of using the algorithms is illustrated 
by several examples. 


32 POLISH MARITIME RESEARCH, No 3/2004 


8. Related and open problems 


Domains of attraction for limit reliability functions of two- 
-state systems are introduced. They are understood as the con- 
ditions which the reliability functions of the particular compo- 
nents of the system have to satisfy in order the system limit 
reliability function to be one of the limit reliability functions 
from the previously established class for this system. Exem- 
plary theorems concerning domains of attraction for limit relia- 
bility functions of homogeneous series systems are presented 
and the application of one of them is illustrated. 

A practically important problem of accuracy of the asymp- 
totic approach to large systems reliability evaluation, which 
concerns the speed of convergence of system reliability sequ- 
ence, is discussed. This problem is illustrated by analysing the 
speed of convergence of the homogeneous series-parallel sys- 
tem reliability sequences, to its limit reliability function. Series - 
-"m out of n" systems and "m out of n"- series systems are 
defined and exemplary theorems on their limit reliability func- 
tions are presented and applied to the reliability evaluation of 
an illumination system and a rope elevator. 

Hierarchical series-parallel and parallel-series systems of 
any order are defined, their reliability functions are determi- 
ned and limit theorems on their reliability functions are ap- 
plied to reliability evaluation of exemplary hierarchical sys- 
tems of 2nd order. 

Applications of the asymptotic approach to large series sys- 
tems reliability improvement are also presented. The chapter 
is completed by showing the possibility of applying the asymp- 
totic approach to the reliability analysis of large systems rela- 
ted to their operation processes. In this scope, the asymptotic 
approach to reliability evaluation for a large port grain trans- 
portation system related to its operation process is performed. 


Miscellanea 


Z IM next tenure 


On 20 May 2004 Marine Technology Unit (acting within 
the Transport Technical Means Section , Transport Commit- 
tee, Polish Academy of Sciences) commenced its next tenu- 
re of activity. ZTM's meeting was held at the Faculty of 
Ocean Engineering and Ship Technology, Gdansk Universi- 
ty of Technology, and chaired, as usually, by Prof. J. Girtler. 


The new tenure was inaugurated 
by presenting two scientific papers : 


# Possible applications of semi-Markovian reliability mo- 
dels to ship power plant operation processes — by J. Rud- 
nicki D.Sc. 

# Probabilistic models in designing the power systems of 
bucket dregers — by D. Bochenski D.Sc. 


Both authors are scientific workers of the Faculty. 


After discussion on the presented topics the second, or- 
ganizational part of the meeting took place. The last-year 
activity of the Unit was summed up and evaluated, the acti- 
vity program for the year 2005 presented, as well as several 
other organizational issues discussed. 


